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The health risks of space radiation are arguably the most serious challenge to space exploration, 
possibly preventing these missions due to safety concerns or increasing their costs to amounts 
beyond what would be acceptable. Radiation in space is substantially different from Earth: high- 
energy ( E ) and charge (Z) particles (HZE) provide the main contribution to the equivalent dose in 
deep space, whereas y rays and low-energy a particles are major contributors on Earth. This 
difference causes a high uncertainty on the estimated radiation health risk (including cancer and 
noncancer effects), and makes protection extremely difficult. In fact, shielding is very difficult in 
space: the very high energy of the cosmic rays and the severe mass constraints in spaceflight 
represent a serious hindrance to effective shielding. Here the physical basis of space radiation 
protection is described, including the most recent achievements in space radiation transport codes 
and shielding approaches. Although deterministic and Monte Carlo transport codes can now 
describe well the interaction of cosmic rays with matter, more accurate double-differential nuclear 
cross sections are needed to improve the codes. Energy deposition in biological molecules and 
related effects should also be developed to achieve accurate risk models for long-term exploratory 
missions. Passive shielding can be effective for solar particle events; however, it is limited for 
galactic cosmic rays (GCR). Active shielding would have to overcome challenging technical hurdles 
to protect against GCR. Thus, improved risk assessment and genetic and biomedical approaches are 
a more likely solution to GCR radiation protection issues. 
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I. INTRODUCTION 

The 20th century will be remembered as a time when 
mankind left planet Earth for the first time. The 21st century 
may well be the “age of maturity” for space exploration. 
Many nations now recognize space as a key driver for educa- 
tional, economic, scientific, and technological advancement 
(ESA, 2007). In addition to the International Space Station 
(ISS), there are four destinations that offer suitable destina- 
tions for human explorers: the Moon, Mars, near Earth 
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FIG. 1 (color). One of the first steps in the 21st century explora- 
tion plans has been the launch of the NASA Lunar Reconnaissance 
Orbiter (LRO) on 18 June 2009. LRO is in low-polar orbit (50 km 
altitude) around the Moon and is designed to generate a compre- 
hensive atlas of the Moon’s features and identifying available 
resources. Measurements of the lunar ionizing radiation environ- 
ment represent one of the main goals of the LRO mission (see 
Sec. II.D.3), and two dedicated instruments are in the payload: 
CRaTER is a cosmic-ray telescope for the detection of charged 
particles with energies >10 MeV, and LEND is a spectrometer to 
measure neutrons at energies <15 MeV. Image: artist’s view of 
LRO courtesy of NASA multimedia gallery. 

objects (NEOs), and the Earth-Moon Lagrangian points. Each 
of these has a unique environment and offers different chal- 
lenges and opportunities. Many robotic lunar missions from 
the USA, Russia, China, Japan, India, and Europe are already 
planned or are under way (Fig. 1). Once the ability to survive 
on the Moon and utilize the alien lunar environment has been 
demonstrated, the way will be open for the first human 
expeditions to Mars. 

Scientific achievements expected from human exploration 
include exobiology and astronomical observatories on the 
Moon or in deep space laboratories. Technological and eco- 
nomic developments are linked to applied microgravity re- 
search, entrepreneurial activities, and space services, such as 
use of the Moon or asteroids as a source of natural resources 
and industrial raw materials, use specific areas of outer space 
as safe havens and operational hubs, and in-space infrastruc- 
tures and services to support space exploration and make it 
sustainable. 

The final goal of space exploration programs is enabling 
life in space through a stable and safe colonization of the 
Solar System. Can this be achieved? Certainly space is not a 
safe environment for humans and enabling life in space 
requires coping with several health problems (White and 
Avemer, 2001). The risks of space travel have been compre- 
hensively summarized in the NASA bioastronautics roadmap 
(NASA, 2005), now updated in the Human Research 
Roadmap 1 , and current gaps in knowledge also identified. 
Risks were rated from 1 (risk of serious health effects, 
and mission could be impossible without mitigation) to 3 

1 http://humanresearchroadmap.nasa.gov/. 


(suspected health consequences with limited impact on the 
mission design). The risks can be summarized into three 
broad categories (NASA, 2009): 

(1) physiological problems caused by microgravity (or 
reduced gravity), 

(2) psychological and medical problems caused by isola- 
tion, and 

(3) acute and late risks caused by exposure to radiation. 

The physiological changes in weightlessness have been ex- 
tensively studied, especially during long-term missions on 
space stations (ISS and, previously, Mir) in low-Earth orbit 
(LEO). Bone loss, kidney stone formation, skeletal muscle 
mass reduction, cardiovascular alterations, impaired sensory- 
motor capabilities, and immune system dysfunctions are 
among the consequences of prolonged permanence in micro- 
gravity. The risks are well characterized and several counter- 
measures are available. None of these risk are given the 
highest rating of 1 (i.e., risk of serious adverse health con- 
sequences, and there is no mitigation strategy validated in 
space or demonstrated on Earth) in the bioastronautics road- 
map (risk category 2 is associated with a serious risk with no 
countermeasure validated in space, and 3 is for known or 
suspected health consequences whose mitigation strategies 
have been validated in space). 

Isolation may lead to serious neurobehavioral problems 
caused by poor psychosocial adaptation (NASA, 2009). 
Several ground platforms are used to study these problems 
and develop countermeasures, such as the Concordia base in 
Antarctica and the Mars500 isolation experiments currently 
under way in Russia. Isolation also brings the problem of 
autonomous medical care (AMC), i.e., the capability to 
handle sickness or accidents in complete isolation. This is 
clearly a risk category 1 for the mission to Mars. 
Countermeasures for AMC risks are mostly technological, 
i.e., rely on the development of portable medical equipment 
and telemedicine. 

Finally, there are the risks related to exposure to space 
radiation (ESA, 2006). Because of the complex nature of the 
space radiation environment (see Sec. II), both acute (i.e., 
short-term risk of radiation sickness) and late (e.g., cancer) 
effects are possible. Acute radiation syndrome can be asso- 
ciated only with intense solar particle events (SPE) with 
crews unable to reach adequate shielding. Late radiation 
morbidity including cancer and other diseases of old age is 
associated with the chronic exposure to galactic cosmic 
radiation (GCR), which is substantially different both quali- 
tatively and quantitatively from the Earth’s radiation natural 
background. Because of the qualitative difference in the 
radiation spectrum (x, y, /3, and a rays on Earth; protons 
and heavy ions in space), the uncertainty in radiation risk 
estimates is very high, especially for carcinogenesis, central 
nervous system (CNS) damage, and late cardiovascular dam- 
age. Early estimates of the uncertainty on space radiation 
cancer mortality risk ranged from 400% to 1500% (NASA, 
1998), with more precise estimates showing uncertainties at 
the 95% confidence level of 4-fold times the point projection 
(Cucinotta, Schimmerling et al., 2001; Cucinotta and 
Durante, 2006). Despite the extensive ground-based experi- 
mental program currently ongoing, this uncertainty is still 
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very high (Schimmerling et al., 1999; Durante and Cucinotta, 
2008). Moreover, countermeasures are not easily available. In 
every radiation protection textbook it is normally stated that 
there are three means to reduce exposure to ionizing radia- 
tion: increasing the distance from the radiation source, re- 
ducing the exposure time, and by shielding. Cosmic radiation 
is isotropic, and therefore increasing the distance is not an 
option in space. Time in space should be increased rather than 
decreased according to the plans of exploration and coloni- 
zation, although reduction of the transit time to the planet, 
where heavy shielding can be more easily achieved, may 
contribute to reducing radiation exposure (Durante and 
Bruno, 2010). Shielding remains as the only feasible counter- 
measure, but it cannot be a full solution for the GCR problem, 
even though it can significantly contribute to risk reduction 
(Wilson et al., 1995). Very heavy shields are impractical on 
spaceships, although small “storm shelters” can be designed 
against intense SPE. Other strategies include the choice of an 
appropriate time of flight, i.e., mission planning and ability to 
predict solar particle events (see Sec. II.B), administration of 
drugs or dietary supplements to reduce the radiation effects, 
and crew selection based on genetic screening. However, it is 
generally acknowledged that the high uncertainty on risk 
estimates and the lack of effective countermeasures make 
cosmic radiation one of the main health concerns for space 
exploration. 


A. Radiation protection quantities 

The physical basis of radiation protection on Earth is well 
established and based on solid scientific evidence and epide- 
miological data, as summarized in the latest report to the 
United Nations (UNSCEAR, 2008) and of the US National 
Academy of Sciences (NRC, 2006). The new recommenda- 
tions of the International Commission on Radiological 
Protection (ICRP) have been published in the recent report 
ICRP-103 (ICRP, 2007). The new ICRP recommendations 
formally replace the commission’s previous (1990) recom- 
mendations (ICRP, 1991). 

As described in Sec. II, the space radiation spectrum is a 
complex mixture of protons, neutrons, and heavy ions with 
energies from a few MeV / nucleon up to TeV/ nucleon. All of 
these components can have substantially different relative 
biological effectiveness (RBE), which is defined as the ratio 
of doses of particles to y rays that produce and equivalent 
biological effect. Therefore, measuring the dose D TR (in gray, 
Gy) averaged over a tissue T by the tissue weighting fraction 
(vv r ) (values shown in Table I) due to radiation R does not 
provide information about the biological response. An ap- 
proximate scheme is to calculate the equivalent dose H T (in 
sievert, Sv) in the organ or tissue T using 

H T = Y^rDt, (1) 


where w R is the so-called radiation weighting factor, whose 
values are based on the review of biological information. The 
revised w R values reported in the ICRP-103 are provided in 
Table I, and the energy-dependent values for neutrons are 
provided in Fig. 2. 


TABLE I. Radiation (w K ) and tissue (w r ) weighting factors, as 
recommended in the latest ICRP-103 report (ICRP, 2007). 


Radiation type 

Wr 

X and y rays 

1 

Electrons and muons 

1 

Protons and charged pions 

2 

a particles and heavy ions 

20 

Neutrons 

2-20 (see Fig. 2) 

Organ 

w T 

Breast, bone marrow, lung, colon, stomach 

0.12 

Gonads 

0.08 

Bladder, liver, esophagus, thyroid 

0.04 

Bone surface, brain, salivary glands, skin 

0.01 

Remainder 

0.12 


Exposure of different organs or tissues is associated with 
different risks of stochastic effects. For example, irradiation 
of blood-forming organs (BFO) will produce a much higher 
risk of leukemia than skin irradiation. If the irradiation is not 
uniform, the different values of H T have to be weighted when 
estimating the risk. Although a whole-body exposure is ex- 
pected in space, doses to different organs can be substantially 
different for solar particle events while very similar for GCRs 
(see Sec. II.D.2). To estimate the risk of all stochastic effects, 
whether the body is exposed uniformly or not uniformly, the 
ICRP recommends the use of a gender and age independent 
effective dose E (also measured in Sv): 


E = 


( 2 ) 


where the indices M and F refer to male and female reference 
phantoms, respectively. Of course, radiation risks depend to a 
large extent on age and gender, making the use of Eq. (2) 
highly inaccurate for risk estimation. There is an inherent 
assumption in the ICRP approach that the worker exposures 
are well below limits making crude approximations 
acceptable. 

In space radiation protection, however, radiation weighting 
factors are not used. As a matter of fact, the space radiation 



FIG. 2. Radiation weighting factor w R for neutrons at different 
energies. The step function is from ICRP60 (1991), and the con- 
tinuous function is the latest ICRP-103 (2007) recommendation. 
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spectrum in different organs is a complicated mixture of 
primary and secondary charged particles with different 
energies. ICRP-103 indeed acknowledges explicitly the 
following: 

“The selection of a single w R value of 20 for all types and 
energies of heavy charged particles seems to be a conserva- 
tive estimate for many particles and is recommended as 
sufficient for general application in radiological protection. 
For applications in space, where these particles contribute 
significantly to the total dose in the human body, a more 
realistic approach may be used.” (ICRP, 2007, para. 124,). 

The approach recommended by the National Council on 
Radiation Protection and Measurements (NCRP), at least for 
LEO (NCRP, 2000), is to use the organ dose equivalent H T 
defined by ICRU (1993) as a mass m average over the tissue: 

Hr = m L dm f Q(L)F T {L)LdL, (3) 

where L is the linear energy transfer (LET = dE/dx, nor- 
mally expressed in keV/ /tm), m is the organ mass, F r is the 
fluence of particles through the organ T, and Q is the quality 
factor (dimensionless), which was adopted in radiation pro- 
tection (ICRP, 1977) prior to the introduction of the weighting 
factors. In radiation protection LET refers to the local energy 
deposition, which is approximated by the stopping power. 
Unlike the weighting factors (Table I), Q is a continuous 
function of LET. The latest recommendations for Q(L) are 
provided by the ICRP-60 (ICRP, 1991) and are given in Fig. 3. 
The integration over the mass of the organ is approximated by 
averaging over a number of points within the organ, for 
example, 33 representative points for the BFO in a reference 
phantom. As discussed in Sec. V, most of the uncertainty in 
Eq. (3) is associated with the quality factor Q. There is also 
clear experimental evidence that the biological effectiveness 
of charged particles depends on both Z and (3, i.e., on the 
track structure, and not on the LET alone (Katz, Scharma, and 
Homayoonfar, 1972; Belli et al., 1993; Cucinotta, Wilson 
et al., 1996; Schardt, Elsasser, and Schulz-Ertner, 2010). 
Nevertheless, Q(L) represents a first approximation of the 



0.1 1 10 100 1000 


LET (keV/qm) 

FIG. 3. Dependence of the quaiity factor Q from LET according to 
the iatest ICRP recommendations. 


changes in RBE in a mixed radiation field such as that in 
space. 

Effective or equivalent doses are ultimately used to esti- 
mate radiation protection limits. The protection for workers 
on Earth relies in the first place on the so-called as low as 
reasonably achievable (ALARA) principle, as formulated by 
ICRP (ICRP, 1977): “All radiation exposures should be kept 
as low as reasonably achievable, economic and social factors 
taking into account.” There is not enough scientific evidence 
to support a lower bound on what constitutes a harmful dose, 
and the issue of low-dose research is always a major compo- 
nent in research programs in the USA and Europe (HLEG, 
2009). In practice ALARA is not a scientific principle, but 
rather a management tool to implement safety factors below 
dose limits using cost-benefit analysis. 

B. Exposure limits 

As yet, terrestrial exposure limits from the ICRP are 
largely based on cancer risk with an added contribution 
from hereditary risks as reflected in the tissue weights of 
Table I. This is due to the evidence that cancer occurs at much 
lower doses than any acute effects. A single coefficient is 
used at all doses assuming a low dose-rate exposure, thus 
following the linear-no-threshold dose-response model, 
which despite many controversies has not been disproved 
(Brenner etal., 2003; Mullenders etal., 2009). The estimates 
of cancer risk attributable to radiation exposure have not 
changed greatly in the past years and continue to be largely 
based on the life span study (LSS) of the survivors following 
the atomic bombings of Hiroshima and Nagasaki. The LSS 
estimates (which correspond to acute exposure) are corrected 
by a dose and dose-rate effectiveness factor (DDREF) for 
solid cancers, which takes into account the reduced radiation 
effectiveness at a low dose rate: ICRP-103 upheld their 
recommendation of the previous DDREF value of 2, while 
the BEIR VII report argues for a lower value of 1.5 which 
would increase solid cancer risk estimates by 33% over the 
older value. Finally, ICRP-103 (ICRP, 2007) still recom- 
mends a risk coefficient for an average population of 
5 X 10 -2 Sv -1 , which was also proposed in ICRP-60 
(ICRP, 1991). The corresponding radiation limits for the 
radiation workers and the general population are summarized 
in Table II. For comparison, typical doses associated with 
medical imaging exposures to ionizing radiation are also 
shown. The average annual effective dose on Earth is about 
3 mSv and is reduced to less than 1 mSv if man-made sources 
and the internal a-particle exposure to Rn daughters are 
excluded. A few areas on Earth display anomalously high 
levels of background radiation, as is the case with Th-rich 
monazite sands, where values 200-400 times higher than the 
world average can be found (UNSCEAR, 2000). It should 
also be noted that approximately 20% of the U.S. population 
receives annual effective doses exceeding the natural 3 mSv 
background (Fazel et al., 2009). Dose rates in space are 
definitely higher than on Earth, although still in the low-dose 
region. The equivalent dose rate in LEO is about 1 mSv/day; 
on Mars it ranges between 100 and 200 mSv/yr (Durante and 
Bruno, 2010), depending on the solar cycle and altitude; 
values around 350 mSv/yr can be reached on the Moon 
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TABLE n. Effective dose limits recommended by ICRP-103 
(ICRP, 2007), and comparison with typical exposures due to medi- 
cal procedures (Fazel et al., 2009). The medical exposures are not 
considered in the computation of the maximum dose. 


Exposure Occupational (mSv) Public (mSv) 

Maximum effective 20 mSv, average over 1 mSv 

dose per year defined periods of 5 yr 

Medical procedure 


Chest x-ray film 

0.02 

(posteroanterior) 


Computed tomography 

7.0 

(CT) chest 


CT abdomen 

8.0 

CT angiography 

15.0 

Mammography 

0.4 

Myocardial perfusion 

15.0 

imaging 


Urography 

3.0 

Nuclear bone imaging 

6.0 


(ESA, 2003). Therefore, astronauts experience a considerably 
higher radiation risk than the population on Earth and occupa- 
tionally exposed workers, particularly considering that the 
latter group receives on average far lower doses than the 
annual limits. NASA has put in place a new approach for 
radiation limits that accounts for uncertainties at the 95% con- 
fidence level (Cucinotta, Hu et al., 2010). This approach 
shows how research to lower uncertainties will be the most 
cost effective method to enable long-term space missions. 

ALARA should be used to optimize radiation protection in 
space missions too. For astronauts, the recommendations of 
NCRP (2000) have been issued only for activities in LEO and 
are based on 3% risk of exposure-induced cancer mortality. If 
the maximum dose for radiation workers in Table II is 
applied, a career lifetime excess risk of approximately 3% 
is indeed reached (assuming a 10 yr career). The risk coef- 
ficients are age and gender dependent, and this was carefully 
considered by NCRP. The corresponding NASA astronaut 
career limits are reported in Table III and compared with 
those from other space agencies. The European, Canadian, 
and Russian space agencies do not adopt age- and gender- 
corrected coefficients, but rather resort to a single career limit 
of 1 Sv. These limits are, however, currently under revision 

TABLE III. Age- and gender-dependent career effective dose 
limits (in Sv) as recommended by different space agencies 
(Cucinotta, Hu et al., 2010; Straube et al., 2010). NASA limits 
are always based on 3% risk of exposure-induced death, and the 
values refer to a 1 yr mission. Different values will be calculated for 
different mission durations (Cucinotta, Hu et al., 2010). 


Space agency 

Gender 

Age at first exposure, (yr) 
30 35 45 55 

NASA (USA) 

Female 

0.47 

0.55 

0.75 

1.1 


Male 

0.62 

0.72 

0.95 

1.5 

JAXA (Japan) 

Female 

0.6 

0.8 

0.9 

1.1 


Male 

0.6 

0.9 

1.0 

1.2 

ESA 


1.0 

1.0 

1.0 

1.0 

FSA (Russia) 


1.0 

1.0 

1.0 

1.0 

CSA (Canada) 


1.0 

1.0 

1.0 

1.0 


(Straube et al., 2010). For exploratory-class missions beyond 
LEO, a description of the major information needed to issue 
recommendations has been instead provided in a dedicated 
NCRP report (NCRP, 2006). 


II. SPACE RADIATION ENVIRONMENT 


The radiation environment in the Solar System is a com- 
plex mixture of particles of solar and galactic origin with a 
broad range of energies. For radiation protection, the prob- 
lems are related to GCRs and particles ejected from the Sun 
during SPEs, as well as secondary radiation produced through 
interaction with planetary surfaces. The solar wind particles, 
even when enhanced due to higher solar activities, do not 
contribute to the radiation burden to astronauts due to their 
relative low energy and hence their absorption in already very 
thin shielding thicknesses. Nevertheless, the solar wind mod- 
ulates the flux of galactic cosmic rays in the energy range 
below about 1 GeV/nucleon, with a fairly regular 11 -yr long 
cycle (see Fig. 4 for current solar activity and predictions for 
the next cycle based on sunspot numbers). During phases of 
higher solar activity the cosmic -ray flux is decreased by a 
factor of 3 to 4 against phases during minimum solar activity. 
However, the probability of SPE is higher during solar maxi- 
mum. The question is then whether a mission should be 
planned in solar minimum (lower chance of acute exposure 
from SPE, but higher GCR risk) or solar maximum. However, 
there is no clear pattern of SPE occurrence with large SPEs 
observed in about 9 of the 1 1 years of any given solar cycle. 

Through the Earth’s magnetic field and an atmospheric 
thickness of about 1 kg/cm 2 thickness, the exposure to cos- 
mic radiation on the Earth’s surface is reduced to a nearly 
zero level. Leaving Earth astronauts are shielded by the 
structure of the spacecraft. For the ISS the interior is shielded 
on average by 20 g/cm 2 , a shielding close to that of the 
Martian atmosphere. In LEO astronauts are still protected 
by the magnetosphere which limits even the exposure to solar 


Solar Cycle 24 Sunspot Number Prediction 



Lov. Prediction (Smoothed) High Prediction (Smoothed) 1 -Sigmo Error 

Smoothed Monthly Values — — Monthly Values 

Updated 2007 Apr 20 NOM/SEC Boulder, CO USA 


FIG. 4 (color). Solar cycle 23 and predictions of the sunspot 
numbers for solar cycle 24 using different models. The next solar 
maximum is expected in late 2001-mid 2012, and the next solar 
minimum around 2019. Data and prediction from the NOAA/NSW 
Space Weather Prediction Center (http://www.swpc.noaa.gov/). 
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energetic particles to exposures far below such causing acute 
radiation effects in man. SPE are, therefore, mostly an issue 
for exploratory-class missions. 

A. Galactic cosmic radiation 

Galactic cosmic radiations originate outside the Solar 
System and impinge isotropically on Earth. Because of their 
high energies (up to 10 20 eV), they most probably originate 
from supernova explosions, neutron stars, pulsars, or other 
sources where high energetic phenomena are involved 
(Cronin, 1999). Detected particles consist of 98% baryons 
and 2% electrons. The baryonic component is composed of 
85% protons (hydrogen nuclei), with the remainder being 
helium (14%) and heavier nuclei (about 1%). Figure 5 
(Cucinotta et al., 2003; Durante and Cucinotta, 2008) using 
the hzetrn code (Wilson et al., 1991; Wilson et al., 1995; 
Tweed, Wilson, and Tripathi, 2004) and the Badhwar-O’ Neill 
GCR model (Badhwar and O’Neill, 1994; O’Neill, 2010) 
shows the relative contribution of the different elements in 
flux, dose, and equivalent dose. The energetic ions heavier 
than helium nuclei have been termed HZE particles (high 
charge Z and high energy E). Although iron nuclei are only 
one-tenth as abundant as carbon or oxygen their contribution 
to the GCR equivalent dose is higher than protons, since the 
dose is proportional to the square of the charge and Q is high 
(Fig. 3). The GCR make up more than 80% of the effective 
doses to crews on ISS, much higher than trapped radiation 
because of their higher penetration power to deep seated 
organs and large quality factors (Cucinotta, Nikjoo, and 
Goodhead, 2000; Cucinotta et al., 2008). 

In addition to the galactic cosmic rays, the so-called 
anomalous component [anomalous cosmic rays (ACRs)] is 
observed. ACRs consist of originally neutral particles coming 
from the interstellar gas which become single ionized by solar 
radiation after entering the heliosphere. These particles are 


FIG. 5 (color). Relative contribution in fluence (circles), dose 
(triangles), and dose equivalent (squares) of different elements in 
the GCR from the HZETRN computer code as reproduced from 
Cucinotta et al. (2003) in Durante and Cucinotta (2008). The 
calculation is an average over 1-year in solar minimum behind 
5 g/cm 2 Al shielding. 


then accelerated in collision regions between fast and slow 
moving streams of the solar wind. They are able to penetrate 
deeper into the magnetic field than fully ionized cosmic rays. 
Their energies are around 20 MeV/nucleon. Therefore, they 
can contribute only to radiation effects behind small shield- 
ing. However, it has to be considered that they lose all then- 
electrons after penetration of a very small amount of shield- 
ing material and thus also deposit energy proportional to the 
square of their charge. 

The propagation of cosmic rays in the interplanetary me- 
dium can be described by a Fokker-Planck equation. If U is 
the cosmic -ray density, E is the particle kinetic energy, V is 
the velocity, and k s is the symmetric part of the diffusion 
tensor, the basic equation is (Parker, 1965) 

— = V(/c s • VU) - V ■ S7U 

+ 5 V# M 1 + ^H-°- (4) 

where m is the proton rest mass. A full numerical solution of 
the equation is given by the deceleration potential cf) (in MV): 

where r B is the radial extent of the heliosphere, k is the 
diffusion coefficient, and V w is the solar wind velocity. The 
deceleration potential is the most important parameter in 
describing the modulation of GCR intensity. An approximate 
solution for the integral fluence j(r, E) at high energies 

300 MeV/nucleon) can be expressed as a function of 
the deceleration potential cf) as 

j(r, E) _ j 0 (r B , E + Ze<t>) 

&-m 2 (E + Zect>) 2 — m 2 ’ 

where y 0 is the local interstellar spectrum. Equation (6) 
basically represents the standard convection-diffusion model 
of the GCR modulation, and although it does not explain the 
radial gradient and the charge dependence, it is widely used to 
predict the fluence rate j(Z, E, t) in any point r in the helio- 
sphere. Since there are no measurements of j 0 , different 
forms for this function have been used, with the constraint 
that the high-energy portion of the spectrum agrees with the 
measurements in LEO. The different GCR models currently 
used basically differ in the choice of this function and the 
solar activity parameter used for prediction. For instance, in 
Nymmik’s model (Nymmik, 1996), also known as the 
Moscow State University model, j 0 is described as a function 
of the particle velocity f3 and rigidity R as 

7o(Z, E)dE = C/B a Ry{^j, (7) 

where C, a, and y are Z-dependent parameters derived from 
fits to experimental data. To obtain the modulated fluence rate 
j n {Z,E,t) at time t during the nth solar cycle, the local 
interstellar spectrum in Eq. (6) must be multiplied by a 
modulation function \P', which is a function of the decelera- 
tion potential, the rigidity, and the Wolf sunspot number. 
Hence, Nymmik’s model is a semiempiric approach which 
relates the solar-cycle variation in the GCR intensity to the 



Charge Number 
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observed time history of the Wolf sunspot number. It is quite 
accurate, with an estimated uncertainty below 15%, and 
therefore commonly used. It is, for example, used in the 
CREME-96 software (Tylka et al., 1997), the most recent 
update of the original cosmic-ray effects on microelectronics 
model developed at the Naval Research Laboratory to predict 
effects of cosmic rays on microelectronics on space. 

There are also other models based on the standard 
diffusion-convection theory of modulation, but different em- 
pirical formulation for the local interstellar spectrum: widely 
used are the CRRES/SPACERAD heavy ion model of the 
environment (CHIME) model (Chen et al., 1994) and the 
Badhwar-O’Neill (or Johnson Space Center) model 
(Badhwar and O’Neill, 1994; O’Neill, 2010). These models 
are based on a fit of the high-energy component of the GCR 
differential flux using a power law: 

it = joP s (e + E 0 )-y, (8) 

where E and E 0 are, respectively, the particle kinetic and rest 
energy and/or nucleon, and y, S, and j 0 are the fitting pa- 
rameters for each charge group. 

The Badhwar-O’Neill model (Badhwar and O’Neill, 1992, 
1994; O’Neill, 2010) was developed upon a self-consistent fit 
of the diffusion equation (4) to a large number of satellite, 
balloon, and more recently data from the advanced compo- 
sition explorer in order to describe solar modulation effects. 
This fit provides the high-energy portion of j 0 , which is then 
modulated using the deceleration potential. The Badhwar- 
O’Neill model has been fit either to the Wolf sunspot numbers 
or to cosmic -ray neutron monitor counting rates as predictors 
of the solar activity. This flux has been measured over long 
periods by different ground-based stations using neutron 
monitors. Despite statistical fluctuations, maxima and min- 
ima of the neutron spectra are clearly anticorrelated to the 
11 -yr solar cycle with a roughly sinusoidal form around an 
average particle flux. 



Energy (MeV/n) 

FIG. 6. GCR particle spectra (H, He, O, and Fe) and their 
modification by solar activity at 1 AU (Badhwar and O’Neill, 
1992; Badhwar, Cucinotta, and O’Neill, 1993) and reproduced in 
NCRP (2000). 


The accuracy of the current models of the GCR is around 
10% for most ions at energies <5 GeV/nucleon, i.e., in the 
region where the solar modulation is most important. Figure 6 
reproduces the differential spectra for H, He, O, and Fe for the 
1977 solar minimum and 1959 solar maximum, the two 
cycles with the strongest modulation of the past 50 years. It 
is noteworthy that the maximum of the spectrum is around 
1 GeV/nucleon for all ions, shifting from lower to higher 
energies and lower fluence in solar maximum. 

B. Solar particle events 

Besides electromagnetic radiation, the Sun continuously 
emits particle radiation, consisting mainly of protons and 
electrons, the so-called solar wind. The intensities of these 
low-energy particles vary from 10 10 to 10 12 particles 
cm _2 s _1 sr _1 with velocities between about 300 and 
800 kms -1 . The related energies are so low (for a proton 
between 100 eV and 3.5 keV) that the particles will be 
stopped within the first few hundred nm of skin. They are 
therefore not of concern for radiation protection. 

However, occasionally the surface of the Sun releases large 
amounts of energy in sudden local outbursts of hard and soft 
x rays and radio waves in a wide frequency band. In these 
SPE large currents and moving magnetic fields in the solar 
corona accelerate solar matter [reviewed in NCRP (2006)]. 
Coronal particles, mostly protons (McGuire, von Rosenvinge, 
and McDonald, 1986) with a small fraction of heavier nuclei 
(Table IV), with energies up to several GeV escape into the 
interplanetary space. They spiral around the interplanetary 
magnetic field lines. Within the ecliptic plane field lines 
expand from the Sun into the interplanetary medium similar 
to the beam of a rotating garden hose. They connect the Earth 
with a certain spot on the western part of the Sun. The number 
and energy spectrum of particles observed in SPE at Earth 
depends on this connection. SPE show an enormous varia- 
bility in particle flux and energy spectra, but the most intense 
events (Fig. 7) have the potential to expose unshielded space 
crews to life threatening doses (Kim et al., 2009a). Since 
high energetic particles are arriving first and followed by 
particles of lower energies the energy spectrum of SPE 


TABLE IV. Elemental composition of SEP (McGuire, von 
Rosenvinge, and McDonald, 1986). The solar energetic particle 
baseline (SEPB) is defined as the ratio of the flux of each element 
compared to oxygen. Protons account for about 98% of the compo- 
sition in flux, but helium ions can give a contribution around 10% in 


Particle 

SEPB 

Particle 

SEPB 

H/O 

3500 ± 500 

Na/O 

0.0083 ± 0.0015 

He/O 

53 ± 5 

Mg/O 

0.183 ±0.010 

Li/O 

<0.001 

Al/O 

0.0115 ±0.018 

Be/O 

<5 X 10" 4 

Si/O 

0.147 ±0.009 

B/O 

<5 X 10" 4 

P/O 

0.0014 ± 0.0006 

C/O 

0.454 ± 0.018 

s/o 

0.0229 ± 0.0025 

N/O 

0.129 ±0.008 

Ar/O 

0.0016 ± 0.0007 

0/0 

1.00±0.031 

Ca/O 

0.0076 ± 0.0016 

F/O 

<5 X 10" 4 

(Ti + Cr)/0 

0.0024 ± 0.0009 

Ne/O 

0.128 ±0.008 

Fe(group)/0 

0.066 ± 0.006 
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FIG. 7 (color). Integral energy spectra of intense solar particle 
events in the 20th century [Kim, Tylka, Atwell, and Cucinotta 
(unpublished)]. 


FIG. 8 (color online). GCR deceleration potential [Eq. (2)] and 
large (proton fluence > 10 8 cm 2 at E > 30 MeV) GLE occur- 
rences as a function of time. From NASA, 2009. 


particles observed at Earth depends on time after the onset of 
the event. Differential energy spectra of solar particle events 
can be represented by a power law (Nymmik, 1996): 

% = f 0 kaE a - l e- kEa , (9) 

dE 

where / is the proton fluence and f 0 , k, and a are three fitting 
parameters. However, other functional forms can be used, 
including a Weibull function or a double power law in the 
rigidity function (Band, 1993; Tylka and Dietrich, 2009). 
Different functional forms to represent SPE spectra will often 
result in different predictions of doses behind shielding. 

A well-connected SPE with high particle fluxes observed at 
Earth is a rare event which is likely to be observed during the 
period of increasing and decreasing maximum solar activity. 
Therefore, major SPEs are observed at Earth as random 
events with a low frequency, typically one per month. They 
last for several hours or days. Events with significant fluxes of 
protons with higher energies can be observed as “ground- 
level events” (GLE) by neutron monitors (Tylka and Dietrich, 
2009). Figure 8 shows the number of GLEs observed over the 
last solar cycles and their correlation with the deceleration 
potential 4> [Eq. (5)], i.e., the solar-cycle phase. Long gaps 
with no events can be seen during solar minimum activity. 
Between the last GLE in cycle 21 and the first one in cycle 22 
there was a 65 month quiet period which was followed by a 
sequence of 1 1 GLEs within 1 yr. 

For radiation protection purposes, it is important to predict 
intensity, energy, and duration of a SPE. In general, energies 
of SPE are lower than GCR, and shielding is therefore a 
possible solution. Storm shelters, i.e., small areas with thick 
shields, are included in spacecrafts and can be built in 
planetary stations. Maximum risk will be during extra- 
vehicular activities (EVA) with minimal shielding, especially 
planetary exploration when astronauts can be very far from 
the Moon (or Mars) base. Since, unlike GCR, intense SPE can 
be life threatening or induce acute radiation syndromes (in- 
cluding nausea and possibly vomiting), they represent not 
only a health problem, but an operational issue of great 


concern. Intense SPE are seldom, but their occurrence is 
not negligible. The large solar flares of August 1972 
(Fig. 7) might have caused prodromal syndromes for an 
unprotected crew on the lunar surface (Hu et al., 2009). 
Only four months before the event, Apollo 16 was on the 
Moon. Although the events in Fig. 7 are generally used to 
construct “worst case scenario” energy-flux spectra for plan- 
ning shielding and protection (Spillantini et al . , 2007), we are 
aware of at least one event in the past 500 years much more 
intense than anyone else recorded: the so-called “solar super- 
storm,” also known as the “Carrington event,” because it was 
recorded by the British astronomer Richard Carrington from 
28 August to 2 September 1859 (Cliver, 2006). Apart from 
the original reports and visual observations of the event, the 
analysis of nitrate concentrations in ice cores proves that 
the protons fluence at energies >30 MeV was close to 
2 X 10 10 cm” 2 at the top of the polar atmosphere, and it 
would have been similar in LEO. Reconstruction of the 
energy spectrum is particularly problematic, but the lack of 
10 Be in the ice core suggests that this was a soft event. Using 
this spectrum, it has been estimated that the skin dose for an 
unprotected crew in EVA would have been about 35 Gy, but a 
normal spacecraft shielding of 5-10 g/cm 2 Al would be 
sufficient to reduce the BFO dose to 0.2-0.7 Gy (Townsend 
et al., 2006). Events of this intensity are very rare, perhaps 
1-2 per millennium, and clearly their effect strongly depends 
not only on the fluence, but on the energy spectrum. A 
Carrington event with a spectrum such as, for instance, the 
September 1989 event (Fig. 7), would lead to much higher 
doses, radiation sickness even for crews inside spacecrafts, 
and catastrophic consequences for the electronics in space. 

Because intense SPE are seldom, of relative short dura- 
tion, and can generally be shielded effectively, any risk 
would be minimized if an adequate forecast system would 
be available. Unfortunately, our current forecasting ability is 
quite limited. 

Several statistical models of SPE occurrences are avail- 
able: the King model (King, 1974), the Jet Propulsion 
Laboratory (JPL) model (Feynman et al., 1993), the 
Emission of Solar Protons (ESP) model (Xapsos et al.. 
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2000), and the NASA SPE Propensity Model (Kim et al., 
2009a, 2009b). 

The King and JPL models are useful for predicting event 
fluences for long-term degradation but do have limitations 
due to the incomplete nature of the data sets upon which they 
were based. For instance, they do not accurately describe the 
complete distribution, allow for the possibility of infinitely 
large events, and lack strong physical and mathematical 
justification. In the ESP model a more accurate approach, 
based on the maximum entropy theory, is used to generate the 
initial distribution of SPE fluences. The basic idea is to 
maximize the entropy S defined for the variable M = log / as 

S= ~ f p{M)\ogp(M)dM. (10) 


Defining the boundary conditions, a system of mathematical 
equations is used to find the solution p(M) that maximizes S, 
by means of the Lagrange multiplier technique. 

The model developed at NASA Johnson Space Center 
(Kim et al., 2009a, 2009b) considered the 370 SPEs that 
were identified during solar cycles 19-23. Fluence data were 
combined over all 5 cycles to estimate an overall probability 
distribution of d> 30 for an average cycle. Between the years 
1561 and 1950, 71 SPEs with d> 30 > 2 X 10 9 protons cm -2 
(McCracken et al., 2001) were also identified from impulsive 
nitrate enhancements in polar ice cores, and these events were 
shown to have a similar probability of occurrence as the large 
events occur during solar cycles 19-23 (Kim et al., 2009a). 

In the Poisson process model, the propensity for SPE 
occurrence at time t is a continuous function of t, known as 
the hazard function. More formally, the hazard function, say 


A(t) = lim 


P{SPE 0 


fence (6 t + h)} 

h 


( 11 ) 


In terms of A (t), E{N(t h f 2 )}, the expected number of events 
in an interval ( t h t 2 ) is equal to m(fj) — m(t 2 ), where m(t) 
known as the “mean value function,” is the cumulative 
hazard: 


( 12 ) 


The functional form best explaining the data was found to be 
proportional to a beta distribution density function offset by a 
quantity A 0 ; i.e.. 


= , _K_ r (p + q) ( t \p- 

4000 4000 r(p)r($) Uoocy 

(0 < t < 4000) 



is used to predict the number of events over a given threshold 
of 30, 60, and 100 MeV proton fluence (Kim et al., 2009b). 

In addition to knowing the expected number of SPEs of a 
given size within a mission period, the actual distribution of 
N a (t 2 ) ~ Njt \ ), the numbers of such events, provides useful 
information for mission planning. Using the fitted Poisson 
model, the estimated cumulative probability distributions of 
N a (t 2 ) ~ N a (t\ ) can be calculated (Kim et al., 2009b): for a 
120-day mission the probability of at least one “large” SPE 
($30 fluence greater than 10 7 protons cm -2 ) is about 0.8 
(80%), with the probability of having at least two large 
SPEs, about 50% and about 20% for more than two SPEs, 
etc. As mission durations reach 540 days, there is almost no 
chance of avoiding an SPE whose event size is greater than 
10 7 protons cm -2 . The probability for multiple numbers of 
SPE occurrences is increased for longer missions and de- 
creased for higher thresholds. 

All these models can be used to estimate the probability of 
large SPE and planning of shielding amounts based on proba- 
bilistic risk assessment concepts, but do not represent a 
“warning system” of practical use. It was recently shown 
(Posner, 2007) that detection of relativistic solar electrons 



for a “typical” nonspecific cycle of 4000 days’ duration, 
where A 0 , K, p, and q are parameters to be estimated. After 
redefining SPE occurrence times relative to 4000-day cycles 
and combining across all 5 cycles, the method of maximum 
likelihood was used to estimate the unknown parameters 
A 0 , K, p, and q (Kim et al., 2009a). From Eq. (10), it 
can be shown that p,, the time of peak hazard, is 
4000 ( p — \ ) / (p + q — 2) days into a cycle. For the observed 
data, /u. was estimated at 1783 days. Using this approach the 
number of events for a given mission length and period 
within the mean cycle can be predicted. A similar analysis 


FIG. 9 (color). This color matrix provides a color code for the 
future proton intensity, 1 hour ahead of time, as predicted by 
relativistic electron measurements. The parameter space is given 
by the current maximum electron increase parameter, going back in 
time for at least 5 min, but up to 60 min, and the current relativistic 
electron intensity. The matrix is derived from the aggregate of all 
1998-2002 relativistic electron observations and their correspond- 
ing 30-50 MeV proton intensities 1 hour later. The color shows the 
average for the proton intensity in each locus. Low statistics limit 
the extent of the matrix to the bottom and upper right. From Posner, 
2007. 
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may enable a prewaming of up to 1 hour of subsequent proton 
events and allow prediction of the integral number of protons. 
The color matrix in Fig. 9 provides a code to predict future 
proton intensity, 1 hour ahead of time, as predicted by 
relativistic electron measurements. Quantities reported on 
the two axes are the current maximum electron increase 
parameter, going back in time 5 to 60 min, and current 
relativistic electron intensity. The matrix is derived from 
the aggregate of all 1998-2002 relativistic electron observa- 
tions and their corresponding 30-50 MeV proton intensities 
1 hour later. The database used contains both “impulsive” 
(flare-accelerated) and “gradual” (coronal-mass-ejection- 
accelerated) events. The color shows the average for the 
proton intensities in each locus. Statistical considerations 
limit the utility of the matrix at the bottom and upper right 
ranges. These findings not only provide evidence of up to a 
1-hour early detection capability, but also may allow astro- 
nauts and mission control to predict if an event is likely to be 
insignificant, which is the most likely outcome. The fore- 
casting potential is limited for poorly connected events with 
very slow increase in proton intensities. However, these 
events are less hazardous and now casting can be more 
useful than forecasting. Based on these observations, it could 
be possible to build on the Moon a solar particle alert 
station, with electron current detectors used to predict the 
arrival and intensity of SPE in about 1 hour (NASA, 2009). 


C. Trapped particle radiation 

Protons and electrons trapped in the Van Allen belts 
(trapped radiation) represent an important contribution to 
radiation exposure in LEO, such as for astronauts on the 
International Space Station (NCRP, 2000), albeit the GCR 
are the dominant component to astronaut organ dose equiv- 
alents (Cucinotta et al., 2000; Cucinotta et al., 2008). 
Trapped radiation originates from the interaction of GCR 
and solar particles with the Earth’s magnetic field and the 
atmosphere. Electrons reach energies of up to 7 MeV and 
protons up to 600 MeV. A small fraction of heavier ions is also 
observed, at energies <50 MeV/nucleon; however, because 
of their limited penetration capacity, they are of no relevance 
for radiation protection. Charged particles with these energies 
moving into a dipole field can never enter into inner areas of 
this field. However, if they are put into this field for any 
reason, they are restricted to certain positions and cannot 
escape. They move in spirals along the geomagnetic field 
lines and are reflected back between the magnetic poles, 
acting as mirrors (Rodger and Clivelrd, 2008). Different 
processes contribute to fill in particles into the radiation belt 
and two main zones of captured particles are observed. The 
inner belt is mainly formed by decaying neutrons, coming 
from the atmosphere in which they are produced in cosmic- 
ray interactions and producing protons and electrons. The 
outer belt consists mainly of trapped solar particles and is 
populated mainly by electrons. During disturbances of the 
magnetosphere by magnetic storms related to SPE where the 
geomagnetic cutoff is usually depressed particles of lower 
energies can penetrate from outside toward the inner regions 
and fill them up. The radiation belts extend over a region from 
200 km to about 75 000 km around the geomagnetic equator. 


The trapped radiation is modulated by the solar cycle: 
Proton intensity decreases with high solar activity, while 
electron intensity increases, and vice versa. Diurnal variations 
by a factor of between 6 and 16 are observed in the outer 
electron belt, and short-term variations due to magnetic 
storms may raise the average flux by 2 or 3 orders of 
magnitude. The center of the inner belt is quite stable, 
especially with respect to protons with a slight westward drift 
observed over time and in some cases transient belts are 
observed (Badhwar et al., 2002). However, at the lower 
edge of the belt, electron and proton intensity may vary by 
up to a factor of 5. For the majority of space missions in LEO, 
protons deliver the dominant contribution to the radiation 
exposure inside space vehicles. Because of their higher en- 
ergies and correspondingly longer range, their total dose 
surpasses that of electrons at mass shielding above about 
0.3 g/cm 2 Al. At lower shielding (e.g., in case of EVA) the 
absorbed dose is dominated by the electron contribution and 
may reach up to 10 mSv/day. 

Of special importance for LEO is the so-called “South 
Atlantic anomaly” (SAA), a region over the coast of Brazil, 
where the radiation belt extends down to altitudes of 200 km. 
This behavior is due to an 11° inclination of the Earth’s 
geomagnetic dipole axis from its axis of rotation toward 
Northern America and a 500 km displacement of the dipole 
center toward the Western Pacific, with corresponding sig- 
nificant reduced field strength values. Almost all radiation 
received in LEO at low inclinations is due to passages 
through the SAA (Fig. 10). At an orbit with 28.5° inclination, 
six orbital rotations per day pass through the anomaly, while 
nine per day do not. Although traversing the anomaly takes 
less than about 15 min and occupies less than 10% of the time 
in orbit, this region accounts for the dominant fraction of total 
exposure on ISS. 

There are several trapped radiation environments available, 
and the Space Environment Information System (SPENVIS) 
platform developed by ESA (Heynderickx et al., 2004) 
implements three proton models and four electron models. 
SPENVIS is a web-based tool providing space engineers and 
scientists with information on the space environment and its 
likely effects on space systems. The NASA AP-8 and AE-8 



FIG. 10 (color). Space radiation field measured onboard the MIR 
space station in 1995 using a tissue-equivalent proportional cham- 
ber. Image courtesy of Gautam Badhwar (NASA). 
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radiation belt models (Vette, 1991) are still the standards for 
engineering applications and are, of course, implemented in 
SPENVIS. Up to now, the NASA codes are the only models 
that completely cover the region of the radiation belts and 
have a wide energy range for both protons (0. 1-4600 MeV) 
and electrons (0.04-7 MeV). The maps are based on data 
from more than 20 satellites from the early 1960s to the mid- 
1970s, but a considerable part of the range of the NASA 
models was achieved by extrapolation. 

D. Radiation measurements in spaceflight 
1. Particle and neutron spectra in low Earth orbit 

The radiation detector systems applied for the monitoring 
of the radiation environment onboard the ISS can be grouped 
into two categories. The first category comprises the “opera- 
tional radiation monitoring devices,” i.e., area monitors and 
personal dosimeters [such as thermoluminescence dosimeters 
(TLD)] used to measure field quantities and organ or tissue 
doses to be used for normalizing radiation transport calcu- 
lations, individual dose assessment and record keeping pur- 
poses, and real-time or near real-time estimates of dose rates 
for purposes of immediate dose management or ALARA. 

The second category consists of instruments selected for 
specific science topics and for limited operation times, typi- 
cally for one or two increments. These topics are a charac- 
terization of the radiation field, such as determination of 
particle flux and energy spectra for each particle type at 
different locations and times in the different modules and 
outside the space station, and influence of shielding material 
composition and thickness. Presently on ISS there are two 
Italian silicon strip detectors, ALTEA (Narici, 2008) and 
ALTEINO (Casolino et al., 2007), a German (DOSTEL) 
and a Bulgarian (LIULIN) silicon telescope, and a Japanese 
Bonner ball neutron detector. As a by-product, these instru- 
ments deliver data which supplement missing information 
from the operational systems and thereby increase the relia- 
bility of such data substantially. Figure 11 provides an ex- 
ample of LET spectrum measured by DOSTEL (Reitz et al., 
2005). A charged particle spectrum measured on ISS by the 
Sileye-3/ALTEINO instruments is shown in Fig. 12. 



FIG. 11. Comparison of measured LET spectra from DOSTEL 
and CR-39 nuclear track detectors. From Reitz et al, 2005. 



Z 


FIG. 12 (color online). Histogram of particle counts showing the 
nuclear identification capabilities of Sileye-3/ ALTEINO from C to 
Fe in the Pirs module of the ISS (Casolino et al., 2007). The peak at 
Z < 6 is an overlap of different light ions. 


Neutron dosimetry is particularly challenging in space 
flights. Neutrons are produced in nuclear interactions be- 
tween trapped protons or GCR with the spacecraft walls 
and the human body. Many neutron detection methods rely 
on the measurement of charged particle produced by neu- 
trons, which is quite difficult because of the many charged 
particles of similar energies produced by trapped protons, 
primary GCR and projectile fragments. They have a very 
wide energy spectrum extending from thermal neutrons to 
more than 1 GeV. TLDs have low sensitivity to neutrons. 
However, chips with 6 LiF (TLD600) or 7 LiF (TLD700) have 
almost the same response to y rays but different response to 
neutrons. Table IV shows the results of neutron measure- 
ments by nuclear emulsions in different shuttle missions 
(Badhwar, Keith, and Cleghorn, 2001). Recoil proton mea- 
surements in nuclear emulsion can measure neutron energy 
up to about 14 MeV, and data in Table IV suggest that these 
neutrons contribute about 7% of the absorbed dose, but 
around 30% of the equivalent dose. The models predict that 
neutron dose rates of nearly 0.5 mSv/day are possible on the 
ISS with lower estimates below 0.1 mSv/day. Considering 
that the charged particle dose rate on the Mir station was 
between 0.4 and 1 mSv/day, it is clear that neutron contri- 
bution can be comparable to the charged particle dose. 
However, an important comphcation for neutron dosimetry 
in space is the possible double counting of secondary radia- 
tion. Since neutrons are secondaries from primary protons 
and heavy ions, and tertiary or higher-order reactions are non- 
negligible, it is not always clear how the results from different 
devices can be added. Of note is that ISS astronaut biodosim- 
etry measurements using chromosomal aberrations which has 
excellent response to neutrons are quite consistent with data 
from tissue equivalent proportional counter (TEPC) and ra- 
diation transport codes, both of which consider charged 
particles and neutron contributions in a consistent manner 
(Cucinotta et al., 2008). 

Measurements of the neutron energy spectra have also 
recently been performed with other instruments, for instance, 
using Bonner balls on ISS (Koshiishi et al., 2007). 
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tained inside the ISS with the neutron measurements inside the other 
spacecraft; the fitted neutron spectrum evaluated from the measure- 
ment inside the STS-28 in 1989, as well as the neutron spectra 
measured inside the MIR in 1991 under 40 g/cm 2 shielding thick- 
ness, from 1990 through 1993 under 30 g/cm 2 shielding thickness, 
and from 1990 through 1992 under 20 g/cm 2 shielding thickness. 
From Koshiishi el al., 2007. 


Comparison of different measurements on the Space Shuttle, 
MIR, and ISS is provided in Fig. 13. Comparison between the 
ISS and the MIR spectra for E > 1 MeV displays that the 
neutron spectrum has higher flux at greater shielding depth, 
demonstrating that a thicker shield produces more high- 
energy neutrons, while spectra are similar at low energy. 
However, the ISS spectrum has a much lower total flux 
than MIR even though the shielding is similar. The discrep- 
ancy is likely to be due to the different instruments used: 
Bonner ball neutron detectors on ISS and nuclear emulsions 
on MIR and the Space Shuttle. The results suggest that the 
data in Table V should be taken with some caution, and 
precise measurements of neutron spectra and fluxes are still 
needed. 

Figure 14 gives a summary of the effective doses received 
by NASA astronauts during various space flights (Cucinotta 
et al., 2008). Roughly 50% of the absorbed dose is caused by 
trapped protons and 50% by GCR (Benton and Benton, 
2001), whereas only 20% of the effective dose is attributable 
to trapped protons (Cucinotta et al., 2008). The effective 
doses vary with altitude and inclination for each flight. The 



FIG. 14 (color). Summary of mission personnel dosimetry from 
all past NASA crews. Effective dose and population average 
biological dose equivalent for astronauts on all NASA space mis- 
sions, including Mercury, Gemini, Apollo, Skylab, Apollo-Soyuz, 
Space Shuttle, NASA-Mir, and ISS missions. Biodosimetry data 
were collected by measurements of chromosomal aberrations in 
peripheral blood lymphocytes before and after the mission, and 
comparison with individual calibration curves measured on the 
ground. From Cucinotta et al., 2008. 

highest values were observed during the high altitude shuttle 
flights at low inclinations with up to 4 mSv/day and during 
the Apollo program with about 3 mSv/day. The radiation 
fields responsible for this exposure are quite different. 

2. Organ doses 

Measurements of the effective dose by Eq. (2) require 
knowledge of the doses in different organs. Excess relative 
risks for radiation-induced cancers from atomic -bomb survi- 
vors are based on bone marrow doses (about 79% of the skin 
dose) for leukemia and colon dose (about 67% of the skin 
dose) for solid cancers (ICRP, 2007). The recommended 
career limits for activities in LEO (Table III) refer to effective 
doses, which must be estimated from personnel dosimetry 
combined with radiation transport codes (NCRP, 2002). Only 
the skin dose is measured in all astronauts by TLD. The use of 
phantom measurements can improve effective dose estimates 
but ultimately cannot represent the movements of individual 
astronauts within a complex spacecraft. 

Organ doses have been calculated by NASA using compu- 
terized anatomical models (CAM) (Billings, Yucker, 
and Heckman, 1973; Atwell, 1994). The CAM man model 


TABLE V. Comparison between dose and dose equivalents for neutrons and charged particles in 
four different STS missions at 28.5° inclination in LEO. Neutron dose was measured by nuclear 
emulsions and charged particle dose by TLD-100 detectors (Badhwar, Keith, and Cleghom, 2001). 


Mission 

Altitude 

(km) 

Neutron 

(AiGy/day) 

Charged particle 
(F-Gy/day) 

Neutron equivalent 
(/iSv/day) 

Charged particle equivalent 
(jU-Sv/day) 

STS-55 

302 

5.9 

57.2 

52.0 

120.1 

STS-57 

470 

25.3 

461.9 

220.0 

859.4 

STS-65 

306 

11.0 

75.2 

95.0 

157.8 

STS-94 

296 

3.7 

101.5 

30.8 

213.9 
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represents a 50th percentile U.S. Air Force male. The com- 
puterized anatomical female model has been obtained by 
removing male organs and incorporating female organs and 
taking into account the smaller stature of the woman. In the 
10-day STS-91 mission (June 1998), organ doses were mea- 
sured using 59 chips (TLD plus the plastic nuclear track 
detector) positioned into a life-size human phantom torso, 
made of tissue-equivalent resin and including a skeleton 
(Badhwar et al., 2002). The effective dose E [Eq. (2); tissue 
weighting factors are provided in Table I] was 4. 1 mSv, which 
is approximately 90% of the skin equivalent dose. Bone 
marrow and colon equivalent doses were less than 80% of 
the skin equivalent dose. A phantom still onboard ISS is 
MATROSHKA (Fig. 15). MATROSHKA (Reitz et al., 
2009) is designed as an anthrophomorphic upper torso 



FIG. 15 (color). MATROSHKA on ISS: The MATROSHKA ex- 
periment facility (upper left panel) was installed on the outside of 
the ISS on 27 February 2004 to measure radiation doses that 
astronauts experience during EVA activities. The phantom and its 
radiation detectors were brought back inside the ISS on 
18 August 2005 as part of Expedition 11 EVA activities: ISS 
Expedition 1 1 Commander Sergei Krikalev (left) and Flight 
Engineer John Phillips (right) dismantle the phantom and remove 
the passive detector systems for return to Earth (upper right panel). 
Since then, MATROSHKA is measuring doses to different organs, 
using active and passive detectors, and different shielding materials 
covering the torso inside of the ISS. Bottom: The median-saggital 
plane of the 3D of MATROSHKA including “skin dose” distribu- 
tion (Reitz et al., 2009). 2D dose distributions for a slice of the 
head, the shoulder region and the lower torso (the whole phantom 
consists of 33 slices, each 25 mm thick loaded with dosimeters). 
These results are based on dose values measured inside the phantom 
and directly on the skin by thermoluminescence dosimeters. This 
represents the dose distribution in an outside exposure. Graphics 
courtesy of Gunther Reitz (DLR). 


phantom, with 33 slices, each 25 mm in thickness. A base 
structure houses the electronic boxes for data acquisition and 
data transmission of the internal temperature and pressure as 
well as the active radiation detector sensors. A carbon fiber 
container with a mass density of ~0.5 g/cm 2 simulates the 
shielding distribution of an astronauts EVA suit. The phantom 
torso is equipped with TLD placed in 354 polyethylene tubes, 
which are accommodated in the 33 slices, enabling determi- 
nation of the absorbed dose and depth-dose distribution at 800 
measurement points in a 2.5 cm x-y-z grid. Combination of 
TLD and nuclear track etch detectors are assembled in poly- 
ethylene boxes (60 X 40 X 25 mm) at selected organ loca- 
tions (eye, lung, stomach, kidney, and intestine) as well as in a 
Nomex-made travel poncho. For the determination of the skin 
dose detectors are sewed into polyethylene stripes directly on 
the skin surface. Seven active radiation detector monitors are 
distributed in the simulated organ regions — 5X silicon de- 
tectors for determination of the heavy ion and neutron com- 
ponent and one silicon telescope (DOSTEL) is located on top 
of the head of the phantom. MATROSHKA uses therefore 
more than 6000 detectors. The MATROSHKA project started 
in January 2004 and has experienced from that time on three 
exposures of detector sets onboard the ISS: one outside 
exposure (MTR-1) and two inside exposures (MTR-2A and 
MTR-2B). In its third exposure phase (MTR-2B) the 
MATROSHKA facility was stored inside the Russian 
Zvezda module from October 2007 until November 2008. 
Measurements for MTR-2 phase B were performed with 
passive and active radiation detectors. Results from MTR-1 
(Fig. 15) show a reduction in the dose in the internal organs, 
due to the body self-shielding, which is though paralleled by 
an increase in the quality factor (Reitz et al., 2009). For 
instance, the skin dose is about 0.94 mGy/day, but it is 
reduced to 0.2 mGy/day in the stomach. However, the qual- 
ity factor increases from 1.7 to 2.5, and therefore equivalent 
dose rates are 1.6 and 0.6 mSv/day, respectively. 

3. Radiation detectors for exploratory-class missions 

There are no measurements of the radiation environment 
on planetary surfaces such as the Earth’s Moon and Mars. 
Recent measurements in the Moon orbit were performed by 
the RADOM spectrometer in the Chandrayaan-1 Indian lunar 
orbiter (Dachev et al. , 2009) and are currently ongoing on the 
NASA LRO mission (Fig. 1) with two dedicated instruments. 
The cosmic ray telescope for the effects of radiation 
(CRaTER) instrument is designed to characterize the global 
lunar radiation environment and its biological impacts 
(Spence et al., 2010). CRaTER is an LET spectrometer, 
consisting of three pairs of thin and thick silicon detectors 
surrounding two pieces of tissue-equivalent plastic and will 
measure the spectrum in a very broad range, with a lower 
limit in energy around 10 MeV. The lunar exploration neutron 
detector (LEND) instrument measures the neutron emission 
flux from the entire lunar surface (Mitrofanov et al., 2010). 
LEND measures thermal, epithermal, and high-energy neu- 
trons up to about 15 MeV, which are produced by spallation 
of GCR nuclei with the lunar regolith. To satisfy the LRO 
mission measurement requirements, LEND includes four 
high pressure 3 He proportional counters and an organic crys- 
tal scintillator within the collimation module, covered with 
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polyethylene and 'OB. The radiation environment in the Mars 
orbit was measured by the Martian radiation environment 
experiment (MARIE) telescope onboard Mars Odyssey, the 
robotic spacecraft launched in April 2001 and still in orbit 
around Mars (Badhwar, 2004). MARIE collected data during 
Odyssey’s cruise and operated in orbit around Mars until the 
large SPE on 28 October 2003 well beyond the mission 
requirement of 90 days of data collection. MARIE has been 
unable to collect data since that time. As the Earth and Mars 
are often out of phase, SPEs at Mars can be connected from 
solar disturbances from different regions of the Sun, or from 
event propagating from distinct magnetic field lines than 
Earth, thus complicating Earth based forecasts of Mars 
events. MARIE observed 14 events from 2001 to 2003 
(Cleghorn et al., 2004), many of which were not observable 
from Earth and others that were, however, revealing different 
characteristics from Mars. After the demise of MARIE in 
October 2003, a neutron detector on the Odyssey spacecraft 
continued to make observation of enhanced count rates at 
Mars (Zeitlin et al, 2010), which allowed further tests of 
predictive models of SPE propagation. 

It can be concluded that radiation monitoring in future 
interplanetary missions will not be substantially different 
from the current situation on ISS. In its recent recommenda- 
tions on future radiation studies for exploration (NCRP, 
2006), NCRP has listed five topics in the area of dosimetry 
and detectors: 

• validation of radiation transport models, 

• characterization of instruments able to provide micro- 
dosimetric spectra, particularly tissue-equivalent pro- 
portional chambers (TEPC), 

• voxel and physical phantoms for the estimates of the 
effective dose, 

• development of spectrometers able to detect indirect 
ionizing radiation under intense fluxes of charged 
particles, 

• and measurements of neutron production in fragmenta- 
tion processes and neutron interaction cross sections. 

III. SPACE RADIATION TRANSPORT 

The description of the passage of high-energy particles 
through matter can be made using Boltzmann-type transport 
equations that treat the atomic and nuclear collisions that alter 
particle energy and types. As an alternative, Monte Carlo 
(MC) computer codes sample from interaction processes for 
individual primaries or their secondary’s to develop histories 
of charged particle passage and energy deposition in materi- 
als. NASA has developed a Boltzmann equation approach for 
HZE nuclei transport denoted as the hzetrn code (Wilson 
et al, 1991), and coupled light ion-HZE transport also solved 
with this method. The hzetrn code with the quantum frag- 
mentation multiple-scattering model (QMSFRG) cross sec- 
tions (Cucinotta, Kim, and Ren, 2006; Cucinotta etal., 2007) 
agrees with spaceflight measurement of organ dose equiva- 
lent to within ±20% (Cucinotta et al., 2000; Cucinotta et al., 
2008). The Boltzmann equation is a deterministic description 
of particle transport. The hzetrn code with broad energy 
beam boundary conditions and the related greentrn code 
with narrow energy beam or monoenergetic boundary 


conditions (Wilson et al., 1991, 1995; Tweed, Wilson, and 
Tripathi, 2004) calculate the average flux of each particle at 
different depths in shielding. These codes have been well 
validated by space measurements in LEO or accelerator 
measurements where good agreement is found. Recent ex- 
tensions of the code have added multigroup methods for 
bidirectional neutron transport and improved convergence 
methods (Clowdsley et al, 2001; Slaba et al, 2010) and 
the addition of pion transport coupled to the GCR sources 
(Blattnig et al, 2004). Earlier versions of the hzetrn code 
(Wilson et al, 1991) considered only a reduced isotopic grid 
of 59 particles. In recent years the hzetrn code considered a 
complete isotopic grid including the isotopic composition of 
the source particles (Cucinotta et al, 2006b), and to consider 
light ion production through nuclear coalescence (Cucinotta 
et al, 2007). Complex spacecraft geometries are handled 
using ray tracing techniques to represent thousands of space- 
craft parts albeit in a bidirectional transport. 

Monte Carlo codes such as GEANT4 (Agostinelli et al, 
2003), phtts (Niita et al, 2006), fluka (Battistoni et al, 
2008), HF.TG- HF.D S (Townsend, Miller, and Gabriel, 2005), 
and mcnpx (James et al., 2006; Aghara etal., 2009) have 
been used for accelerator and atmospheric radiation studies 
and more recently for space applications. The MC approach 
allows for very detailed considerations on detector designs 
and background events in high-energy experiments to be 
considered. In principle, more complicated geometries can 
be considered in the MC approach compared to the hzetrn 
code because angular scattering of particles is considered; 
however, in practice simplified geometries have been used by 
MC codes to represent complex spacecraft geometries, such 
as cylinders of appropriate dimensions and composition. 
Many MC transport codes include modules to treat low- 
energy neutron and photon transport that were originally 
developed as a reactor design such as the MORSE code. 
There are quite a few nuclear models utilized within the 
different MC codes, including commonalities with nuclear 
models used in the hzetrn codes, which we discuss below. 
What often is overlooked is that the MC codes have the 
ability to evaluate stochastic events (Cucinotta et al., 
2011); however, in most space radiation applications, the 
MC codes have reported only on the calculation of average 
quantities such as fluence, dose, and dose equivalent, and very 
few analyses of the importance of fluctuations are described. 

Space radiation problems are not likely handled with a 
one-size fits all approach and the specific application drives 
the method to be used. One aspect of codes that can be 
standardized is the nuclear interaction models including de- 
veloping methods of cross comparisons of distinct transport 
codes. Computational speeds are optimized if data libraries 
for cross sections are stored and called from data files in 
programs, rather than including the actual calculation proce- 
dure for cross sections within the transport code in order to 
reduce CPU times. Other considerations are the angular 
dependence of the reactions. For example, heavy ions above 
100 MeV/nucleon travel predominantly in the forward di- 
rection and there is only a minor need to evaluate angular 
deflections for the omnidirectional fields in space. For lower 
energy heavy ions, atomic collisions dominate which are well 
described by multiple-scattering and straggling distributions. 
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On the other hand, neutrons are created by several physical 
processes that have distinct angular production characteris- 
tics. Light charged particles will have behaviors intermediate 
to neutrons and heavy ions, and the roles of angular deflec- 
tions for space applications require further study. 

A. Particle and heavy ion transport calculations 
1. Deterministic codes 

The relevant transport equations are derived on the basis of 
conservation principles (Wilson et al., 2001) for the flux 
density cf)j{x, Cl, E) of type j particles as 

Cl • S7(f>j{x, Cl, E) 

= xf cr M n ’ n '’ E ’ E ')4>k^ o', E')dE'dCl' 

- o-j{E)(/)j(x, Cl, E), (14) 

where crj(E) and o-j k (Cl, Cl', E, E') are the media macroscopic 
cross sections. The crj k (Cl, Cl', E, E') represent all those pro- 
cesses by which type k particles moving in direction Cl' with 
energy E' produce a type j particle in direction Cl with energy 
E. The flux density 4>j(x, Cl, E) is the main physical quantity 
used to determine the physical or biological response by 
folding it with an appropriate response function for the 
physical or biological system under study. 

There may be several reactions which produce a particular 
product, and the appropriate cross sections for Eq. (14) are the 
inclusive ones. The total cross section o-j{E) with the medium 
for each particle type of energy E may be expanded as 

o-j(E) = af(E) + af(E) + (r){E), (15) 

where the first term refers to collision with atomic electrons, 
the second term is for elastic nuclear scattering, and the third 
term describes nuclear reactions. The microscopic cross sec- 
tions and average energy transfer are ordered as follows: 

af{E) ~ 10” 16 cm 2 for 8E at ~ 10 2 eV, 

of(E) ~ 10” 19 cm 2 for 8E el ~ 10 6 eV, (16) 

a)(E) ~ 10” 24 cm 2 for 8E t ~ 10 s eV. 

This ordering allows flexibility in expanding solutions to the 
Boltzmann equation as a sequence of physical perturbative 
approximations (Wilson et al., 2001). Many atomic colli- 
sions (~ 10 6 ) occur in a cm of ordinary matter, whereas ~ 10 3 
nuclear Coulomb elastic collisions occur per cm. In contrast, 
nuclear reactions are separated by a fraction to many cm 
depending on energy and particle type. For neutrons, 
crf(E) ~ 0 and the nuclear elastic process appears as the 
first-order perturbation. Mean free paths for elastic scattering 
of neutrons may become quite small, especially at low en- 
ergies in the resonance region (ICRU, 2000). 

The solution of Eq. (14) involves hundreds of multidimen- 
sional integro-differential equations which are coupled to- 
gether by thousands of cross terms and must be solved self- 
consistently subject to boundary conditions ultimately related 
to the external environment and the geometry of the astro- 
naut’s body and/or a complex vehicle. A series of approxi- 
mate solutions can be studied and indicates a high level of 


accuracy for most applications (Wilson et al., 2001; Tweed, 
Wilson, and Tripathi, 2004). The mean energy loss can be 
introduced in a continuous slowing down approximation 
(CSDA). The highly directional Coulomb cross section for 
charged ions and nuclear elastic scattering for neutrons gen- 
erally dominate the second perturbation term. The angular 
dispersion and its effects on lateral beam spread and range 
straggling are important corrections in comparing to labora- 
tory measurements. The nuclear elastic scattering is espe- 
cially important to neutron fields and has been treated in the 
past using Monte Carlo methods or multigroup methods 
(Hughes, Prael, and Little, 1997). The third perturbation 
term consists of complex energy and angle functions. 
Results from Monte Carlo codes (Alsmiller et al., 1965) 
provided the basis for the generation of analytical techniques 
and the simplification of boundary conditions used in space 
shield code development (Wilson et al., 1991). 

For transport of the GCR heavy ions (A > 4), the use of the 
CSDA in a marching procedure has been implemented in 
NASA’s hzetrn code (Wilson et al., 1991). For GCR 
transport, an isotopic grid of more than 150 ions is needed 
to represent all the possible fragments produced with 190 ion 
uses in the hzetrn code (Cucinotta, Kim, and Ren, 2006). 
For monoenergetic beam transport, a Gaussian model of the 
fragment single differential cross section has been considered 
along with energy straggling using a Green’s function ap- 
proach (Tweed, Wilson, and Tripathi, 2004). For light particle 
transport, the broad redistribution in energy of the ions in 
collisions is considered in the marching procedure (Wilson 
et al., 1991). Angular effects for neutron transport are con- 
sidered using a multigroup or Monte Carlo transport models 
(Hughes, Prael, and Little, 1997). 

2. Transport coefficients 

The transport coefficients describe the atomic, molecular, 
and nuclear processes by which the particle fields are modi- 
fied by the presence of a material medium (Wilson et al., 
2001). As such, basic atomic and nuclear theories provide the 
input to the transport code databases. The first-order physical 
perturbation on the right-hand side of Eq. (14) is the atomic 
and/or molecular cross sections as noted in Eq. (16) for which 
those terms in Eq. (14) are expanded about the energy mo- 
ments as 

S n {E) = (17) 

where e, is based on the electronic excitation energy, and 
< Ti{E ) is the total atomic and/or molecular cross section for 
delivering s, energy to the orbital electrons (including dis- 
crete and continuum levels). The first moment (n = 1) is the 
usual stopping power, and the usual CSDA is achieved by 
neglecting the higher-energy moments. 

In Eq. (17) specification of e, and (Ji{E) requires complete 
knowledge of the atomic and molecular wave functions. 
Stopping power databases are derived semiempirically as 
the Bethe reduction of Eq. (17) in terms of mean excitation 
energies and shell corrections (Fano, 1963; Wilson et al., 
1991). The stopping power S is adequately described by the 
Bethe-Bloch formula for most ion energies (Bichsel, 1992): 
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+ Z P L X {$) + Z 2 p L 2 (/3) + Ls(J3 ) }, (18) 


where e is the electronic charge, N T is the density of target 
atoms, m is the mass of the electron, c is the speed of light, 
/3 = v/c , and I is the mean excitation energy. In Eq. (18), the 
various terms are the shell correction C(/3), Barkas correction 
L|(/3), Bloch term L 2 (/3), and Mott and density corrections 
L|(/3). The range of the ion is evaluated from the stopping 
power as 


’) = P 

Jo 


•E dE' 

fey 


(19) 


The second energy moment is related to energy or range 
straggling and provides corrections to the ion slowing down 
spectrum (Fano, 1963). Straggling has been well studied for 
ion beam applications with leading order and most higher- 
order correction terms well understood. For broad beams 
conditions of GCR transport straggling effects are negligible; 
however, they are important for laboratory studies with 
monoenergetic beams and for understanding radiation detec- 
tor responses. The next physical perturbation term is the 
Coulomb scattering by the atomic nucleus and is typically 
represented by Rutherford scattering modified by screening 
of the nuclear charge by the orbital electrons using the 
Thomas-Fermi distribution for the atomic orbits. The total 
nuclear Coulomb cross section found by integrating over the 
scattering directions is related to the radiation length. The 
differential cross section is highly peaked in the forward 
direction, and only after many scatterings is significant 
beam divergence seen. Numerical solutions to the Coulomb 
multiple-scattering problem have been investigated for many 
years (Fermi, 1940) and accurately describe experimental 
data with HZE (Wong et al., 1990) or proton beams 
(Carlsson and Rosander, 1973). 

For determining the particle spectra, numerical techniques 
have been developed which allow for computationally effi- 
cient and accurate computer codes in the straightforward 
approximation (Wilson et al., 1995). In this model the heavy 
ion flux (f) j(E, x) of an ion j with mass number Aj, charge 
number Zj, and energy E (in units of MeV/nucleon) at 
shielding depth x (in units of g/cm 2 ) is determined by the 
partial differential equation (Wilson et al., 1991) 

= ^a hk (E)(l> k (E,x), (20) 


Xj(v ptj) = SjiE^jix.E) (22) 

is introduced leading to the transport equation 

[2-^- + trj'lxiVp €j . ) = ^.Vj.k — xiVb ik), (23) 

L drjj J J k v k 

where S j(E) is the stopping power of ion j, Vj = Z 2 / A a j is 
the total absorption cross section, and is the fragmenta- 
tion cross section to produce isotope j from projectile isotope 
k. The solution of the transport equation for light ions and 
neutrons is distinct from Eq. (23) because of the broad 
redistribution of energy in collision events as described by 
Wilson et al. (1991). 

3. Monte Carlo codes 

There are several MC codes that have been developed for 
accelerator applications in detectors and therapy with high- 
energy charged particles (Durante and Bruno, 2010; Schardt, 
Elsasser, and Schulz-Ertner, 2010), and more recently medi- 
cal and space applications. These codes often differ by 
nuclear interaction models used and scoring and geometry 
approaches. The fluka (fluctuating kascade) code transports 
various charged particles including heavy ions, neutrons, 
electrons, and photons up to 10 4 TeV for all particles and 
down to thermal energies. Particle interactions are scored on a 
3D mesh for event-by-event histories or performing user 
specified averages, hetc-heds is a NASA developed exten- 
sion for the former NASA human exploration and develop- 
ment of space (HEDS) efforts of the HETC (high-energy 
transport code) originally developed at Oak Ridge National 
Laboratory, phits (particle and heavy ion transport code 
system) is a MC code developed in Japan that treats with 
all particles up to 200 GeV/nucleon, and incorporates the 
MCNP code for low-energy neutron and photon transport. 
GEANT4 (geometry and tracking code) is an MC code devel- 
oped at CERN for accelerator and detector applications. 

The MC approach allows more flexibility compared to 
deterministic codes for introducing new applications of in- 
terest, such as particle transport in magnetic fields, and three- 
dimensional aspects that would be prohibitive to solve in an 
analytic approach. Nuclear interaction considerations in the 
various codes are described below. There are distinctions in 
how different MC codes treat complex geometries compared 
to deterministic codes, and much longer CPU times are 
needed for MC codes making them prohibitively time con- 
suming for spacecraft shielding design studies. Many of the 
MC codes provide training and support functions through 
international collaborations. 


where (Xj{E) is the energy-dependent absorption cross section 
(cm -1 ) and <Tj k (E) is the fragmentation cross section for 
producing an ion j from k. The solution to Eq. (20) in the 
hzetrn code is found using the methods of characteristics 
where the coordinate transformation from depth and particle 
range Rj(E) are defined 

Vj = x- Rj(E), Zj = x + Rj(E), (21) 

and the scaled flux 


B. Nuclear interaction cross sections 

The types of cross sections required for transport involve 
total yields and multiplicities and inclusive secondary 
energy spectra for one-dimensional transport or inclusive 
double-differential cross sections in angle and energy for 
three-dimensional transport. For Monte Carlo simulations, 
exclusive cross sections may be needed for computer algo- 
rithms, an enormous construction task when one considers the 
large number of projectile-target combinations, secondary 
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multiplicities, etc. needed to transport all GCR particles and 
energies through spacecraft and tissues. Fortunately, physical 
considerations lead to great simplifications allowing inclusive 
cross sections to be appropriate for most applications. Low- 
energy evaporation products including heavy ion target frag- 
ments are high LET events. Knockout products from proton 
or neutron reactions and projectile fragments from GCR 
nuclei are typically of low to moderate LET; however, their 
large ranges lead to radiation buildup through further 
reactions. 


1. Nuclear fragmentation models 

Three types of nuclear fragmentation models have been 
developed to consider heavy ion fragmentation cross sections. 
The first type uses the multiple-scattering series approach of 
Watson (1953) and high-energy approximations in a quantum 
multiple-scattering approach. The second approach is the 
intranuclear cascade (INC) models, originally developed by 
Metropolis, Bertini, and others at Oak Ridge, and in recent 
years the Liege INC model developed by Cugnon, and one 
based on the patron model for the ultrarelativistic case called 
JAM for jet-AA microscopic transport model. A third ap- 
proach is the molecular dynamics model (QMD), which is a 
semiclassical model using Gaussian wave packets and 
Newtonian equations, and the related Japanese quantum mo- 
lecular dynamics model used by the phits code. Each of these 
models relies on the two-step picture, called cascade and 
evaporation, or abrasion and ablation to describe the interac- 
tion. The description of nuclear reactions through abrasion 
(particle removal during ion-ion interaction) and ablation 
(nuclear deexcitation after the abrasion step) is illustrated in 
Fig. 16, which shows the roles of projectile overlap, fireball 
formation in central regions, and the decay of the prefragment 
spectators. Peripheral collisions lead to small mass removal, 
while central collisions can lead to the total destruction of the 
two nuclei. Ablation applies best to peripheral collisions 
where the remaining nuclei after the collision called the 
projectile or large prefragment are left in a state of excitation 
and will decay to the ground state by statistical emission of 
light particles and y rays. 

The models developed at NASA rely on the quantum 
multiple-scattering theories (QMST), which are solved to 
calculate the probability of abrading n nucleons, and to 
evaluate the excitation spectra of the prefragments. The 
multiple-scattering series of Watson (1953) for proton- 
nucleus scattering was extended to the nucleus-nucleus case 
by Wilson (1974) and forms the basis for considering the 
elastic channel for determination of the total absorption cross 
section, and the development of models for the various 
inelastic channels including nuclear abrasion where the 



Target Target-fragment 



Glauber model (Czyz and Maximon, 1969; Glauber, 2006) 
is useful (Hufner, Schafer, and Schumann, 1975). The equa- 
tions of motion for nuclear scattering are expressed in terms 
of the transition operator which represents an infinite series 
for the multiple scattering of the constituents of the projectile 
and target nucleon. The strong nature of the nuclear force 
requires a nonperturbative solution to the scattering problem. 
In relativistic field theory, the non-Abelian nature of the 
strong force has precluded a formulation of the transition 
matrix for nuclear scattering using the Lagrangian of quan- 
tum chromodynamics (QCD). A relativistically covariant 
formulation of the problem has been developed by Maung 
and co-workers using meson exchange theory (Maung, 
Norbury, and Kahana, 1996). The basic approach, in both 
relativistic and nonrelativistic multiple-scattering theories, is 
to resum the multiple-scattering series, which is expressed in 
terns of the irreducible and reducible exchange diagrams in 
the relativistic multiple-scattering theory (RMST) or the 
nuclear potential in the NRMST, in terms of the transition 
matrix for the constituents of the projectile and target nuclei. 
This avoids having to deal directly with the highly singular 
behavior of the nuclear potential at short distances, and 
instead the constituent transition matrix is used, which is 
often known from experimental determinations. 

The NRMST is obtained by approximating the full Green’s 
function by the leading order tern corresponding to one- 
meson exchange diagrams and using a nonrelativistic reduc- 
tion of the three-dimensional Green’s function. The potential 
term is then the sum of the interactions of the constituents 
(Wilson, 1974) 

V= f.l \ T ‘i (24) 


and the nonrelativistic Green’s function is given by 

g NR = (E — H P — H t )~\ (25) 

where H P and H r are the projectile and target internal 
Hamiltonians, respectively. At high energies the impulse 
approximation is invoked, which assumes that the relative 
kinetic energy of the constituents is much larger than the 
binding energies such that the propagator is given by 

g 0 = (E - T P - T t )~' (26) 

and the constituent interactions are replaced by the free 
interactions which are truly of the two-body form. For 
high-energy reactions, the scattering is often confined to the 
forward direction. Here the eikonal approximation is useful 
for reducing the scattering problem to a closed form expres- 
sion, which reduces the MST to the form of the Glauber 
model (Czyz and Maximon, 1969; Wilson, 1974; Glauber, 
2006), where cancellation of reflection terms occurs, where 
nucleons from the projectile and target can rescatter. 

2. Elastic and inelastic channels 


Abrasion Ablation 

FIG. 16 (color). Illustration of the abrasion-ablation model. 


Knowledge of the elastic amplitude alone is used to de- 
termine the total and absorption cross section. The absorption 
cross section is the key parameter used by MC transport codes 
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to determine if an interaction has taken place in developing 
Monte Carlo histories. The total (TOT) cross section is found 
from the elastic amplitude f{q) by using the optical theorem: 
477 

fr TO T = — Im/(q = 0). (27) 

The absorption cross section is then found as the difference 
between the total and elastic cross sections 

°ABS = a TOT - 0EL> (28) 

where cr EL is the total elastic cross section, and (r A(5S is the 
total absorption cross section. The absorption cross sections 
are accurately represented by energy-dependent variants 
(Townsend and Wilson, 1986) of the Bradt-Peters equation 
(Bradt and Peters, 1950) 

^abs = 7Trlc(EUf + 4 /3 - c 2 {E)f, (29) 

where r 0 , c x (E), and c 2 (E) are parameters. Figure 17 shows 
results from experiments for several projectiles on aluminum 
targets versus the projectile energy. Absorption cross sections 
have been well studied both experimentally and theoretically 
and are known with a few percent accuracy (Tripathi, Wilson, 
and Cucinotta, 2001). Tripathi, Wilson, and Cucinotta (2001) 



Energy (A MeV) 



Energy (A MeV) 


FIG. 17. Comparison of experimental data for p-Al (upper panel) 
or C-Al (lower panel) absorption cross sections to the theoretical 
model of Tripathi, Wilson, and Cucinotta (1998). 


made extension of Eq. (29) to represent most available data 
and this model is used by GEANT4, phits, and hzetrn codes. 

3. Light ion interactions 

For proton or neutron induced knockout of nucleons or 
light clusters and the fragmentation of light nuclei, such as 
cosmic-ray helium, the multiple-scattering series is written in 
a three-body form in order to include the effects of final state 
interactions (FSI) in the formalism, which was developed for 
the hzetrn code (Cucinotta, Townsend, and Wilson, 1992; 
Cucinotta and Dubey, 1994; Cucinotta et al., 1996). An 
alternate approach would be to consider a distorted wave 
form for knocked-out particles. The three-body approach of 
the Faddeev formalism (Faddeev, 1956) allows for a conver- 
gent series if the transition matrix for the light particles 
incident on nuclei is known. The leading order terms to the 
knockout series are represented by overlap functions for the 
virtual decay of the nucleus into the participant and spectator, 
and the quasielastic scattering of the participants. For heavy 
ions the overlap functions are represented by single particle 
wave functions. 

The transition matrix for the reaction P + T —< > a + b + X 
can be written as a three-body problem of a — T, b — T, and 
a — b interactions when rearrangement channels are ne- 
glected and with the understanding that all target final and 
intermediate states must be summed. Using the Faddeev 
method leads to a multiple-scattering series generated by a 
coupled set of integral equations (Cucinotta, Townsend, and 
Wilson, 1992) 

Tl = 1'aT + f bT + 7' ab , (30) 

where 'T a r, 1" br> and I" ab are the cluster amplitudes which 
are the transition operators for aT, bT, and ab scattering, 
respectively, in the projectile-target Hilbert space and G 0 is 
the Green’s function in the impulse approximation. The 
leading order corrections to the pole approximation are found 
by truncating Eq. (30) as 

T" = (1 + '^ab G o){T'aT + T bT + T' gtGqT' b T 

+ 'l' bT G 0 'l' aT ) (31) 

and replacing 'T al and 'T br by their on-shell values. 
Equation (31) allows for all orders of multiple scattering, 
however, assuming the dominance of the ab cluster in the 
projectile and that ab FSI occur only after interactions with 
the target. Comparison of the light ion breakup model to 
experiment for 3 H production a- n C reactions is shown in 
Fig. 18. The MST has been applied to a -particle knockout, as 
well as to the knockout of nucleons from target nuclei by 
incident nucleons. 

Evaluation of the inclusive cross sections for the a-particle 
breakup or light particle knockout from nucleon induced 
reactions involves the quasielastic scattering of the fragments 
(Cucinotta and Dubey, 1994; Cucinotta, Townsend et al., 
1996). For different reactions, the identity of secondaries may 
be the same (e.g., for p and n production) and quasielastic 
scattering of the incident nucleon ( p or n) may overlap with 
the knockout distribution. These individual contributions are 
written as 
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Proton exchange 



p (GeV/c) 

FIG. 18. Transverse momentum distributions for 4 He fragmenta- 
tion on C target into tritons at 1.9 GeV/nucleon. From Cucinotta, 
Townsend, and Wilson, 1992. 


(dv\ _/dc rt + + /*T\ 

\dp/ Pi n Wp/K-O. \dpjQE \dp/ EVAP 

where we also included a contribution from the decay of 
highly excited target recoils. Coalescence of light nuclei 
also should be considered for AA collisions. The quasielastic 
terms may have a contribution from charge-exchange or 
nucleon resonances. 

In the abrasion-ablation models (Hufner, Schafer, and 
Schurmann, 1975; Townsend and Wilson, 1986; Wilson, 
Tripathi et al., 1995; Cucinotta et al., 1998; Cucinotta, 
Wilson et al. , 2006) used to describe heavy ion fragmentation 
cross sections the abrasion and ablation processes are largely 
treated as independent processess only linked by the excita- 
tion energies of the prefragments at the end of the abrasion 
stage. The individual steps of abrasion and ablation can be 
described in both semiclassical and quantum mechanical 
approaches. Nuclear database development has focused 
largely on QMST or Monte Carlo approaches to nuclear 
reactions using an intranuclear cascade model. 

For high energies and increasing charge of the projectile 
and targets, electromagnetic dissociation can lead to the 
knockout of neutrons, protons, and light nuclear clusters. 
The virtual photons created by the electromagnetic field 
of the two nuclei are the cause of this phenomenon. 
Contributions to fragmentation cross sections can be eval- 
uated using the Weizsacker- Williams method of virtual 
quanta where the virtual photon spectra of the target nuclei 
is folded with the photonuclear reaction cross section for the 
projectile emittal (Norbury and Townsend, 1990; Wilson 
et al., 1991). 

4. Intranuclear cascade models 

Monte Carlo simulation techniques can be used to describe 
nuclear multiple scattering (Ferrari and Sala, 1996; Cugnon, 


Volant, and Vuillier, 1997) in the so-called intranuclear 
cascade model (Metropolis et al., 1958). The Monte Carlo 
approach to reaction theories relies on the phase space con- 
siderations and two-body (NN) cross sections and production 
cross sections for various mesons and delta resonances. The 
MST is treated in an algorithmic manner by following each 
nucleon by sampling over the possible energy and momentum 
transfer allowed by the two-body cross sections, often includ- 
ing Pauli blocking and other nuclear medium effects. Several 
approaches are used in these codes (Cugnon, Volant, and 
Vuillier, 1997; Ballarini, Battistoni, and Brugger, 2007) to 
couple to the intranuclear cascade including the preequili- 
brium models and the nuclear evaporation models One ad- 
vantage to the MC approach is that all orders of scattering can 
be followed in the algorithm. However, many of the quantum 
aspects of the problem must be ignored including shell 
structure, the nuclear surface, and interference effects. At 
high energies, the dual patron model implemented in the 
JAM code is used by several of the MC transport codes. 
The fluka code uses a multifaceted approach where at 
high energies the dual patron model is assumed and supple- 
mented by a Glauber-Gribov cascade approach. In the inter- 
mediate energy range a generalize INC model is used and 
finally the preequilibrium approach to nuclear thermalization 
(PEANUT) model is used (Ballarini, Battistoni, and Brugger, 
2007). 

5. QMD 

In the QMD model (Neise et al., 1990) each nucleon is 
represented by a Gaussian wave packet and the total wave 
function is a single particle product of these functions. The 
one-body distribution function is found by the Wigner trans- 
formation of the wave function 

f(r, p) = N exp|^ - ^ ^ -2L(p- P t ) 2 J, (33) 

where L is a parameter representing the spatial spread of the 
wave packet, R and P correspond to the center-of-mass 
position and momentum of the wave packet, respectively, 
and N is the normalization constant. The value of L is taken 
as 2 fm 2 to best represent stable nuclear ground states (Niita 
et al., 2006). The equations of motion for the time develop- 
ment of the center-of-mass coordinates obey Newtonian 
equations. The nuclear potential use in the QMD is usually 
the Skyrme-type NN interaction with a Coulomb potential. 
At the end of the dynamical stage of the reaction, deexcitation 
or ablation of the project and target nuclei occurs using the 
GEM code in the PHITS version of QMD. An arbitrary parame- 
ter in the simulations is the switching time from the dynami- 
cal of the QMD stage to the deexcitation which is usually 
taken as 100 to 150 fm/s (Niita et al., 2006). 

6. Fragment momentum distribution 

The momentum distribution of projectile fragments can be 
described as a Gaussian distribution in the projectile rest 
frame with a small downshift in the average momentum 
from the projectiles velocity. The longitudinal momentum 
width is well described by (Goldhaber, 1974) 
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a L — °0 


p 


Mp - »)] 

A P - 1 J 


1/2 


(34) 


where n is the number of nucleons removed from the projec- 
tile, and cr 0 is related to the Fermi momentum of the projec- 
tile p F by <r 0 = Pf/J 5. The functional form of the model of 
Eq. (34) has been discovered in several distinct models 
(Hufner, 1985); however, it likely arises to the large number 
of intermediate states in abrasion and ablation, such that the 
central limit theorem leads to a Gaussian form. The transverse 
width is approximately the same as the longitudinal for 
heavier fragments. A small momentum downshift also occurs 
and is dependent on the fragment mass. Transformation of the 
Gaussian distribution to the laboratory rest frame at the high 
energies of interest for space radiation transport problems 
reveals a narrow angular distribution for projectile heavy ion 
fragments that are strictly forward peaked in a narrow cone 
(< 5 deg), i.e., the physical reason for the success of the 
straightforward approximation. For lighter fragments the lon- 
gitudinal and traverse widths diverge and the Gaussian model 
breaks down. This is due both to their smaller mass and to the 
multiple sources for light particle production including pro- 
jectile abrasion, projectile ablation, target abrasion, and target 
ablation, as well as a possible intermediate source due to the 
formation of an intermediate rapidity fireball in central 
collisions. 


C. Nuclear fragmentation measurements 

Older studies used chemical decay methods to extract 
fragmentation cross sections from target decays; however, 
these data are considered to be less accurate than the mea- 
surements of the cross section directly from projectile inter- 
actions with target atoms. Nuclear plastic track detectors 
became commonly used for fragmentation cross-section mea- 
surements through the introduction of CR-39 solid state 
nuclear track detectors (C 12 H 18 0 7 ), which allowed lower 
charge ions to be measured increasing the accuracy of frag- 
mentation measurements and this approach is still used today. 
Iancu, Flesch, and Heinrich (2005) used CR-39 detectors of 
0.6 mm thickness to determine trajectories and charges of 
nuclear fragments. Etching in 6 N NaOH at 60° C for 59 h 
helped to identify etch pits and latent tracks which are related 
to the energy loss of the ions. Scanning systems on computed 
controlled microscopes and image analysis software have 
been developed to improve data analysis. Of interest is a 
similar approach that can be used to understand charge and 
FET spectra from astronauts returning from space when 
CR-39 is made part of crew personnel dosimetry. 

The earliest studies of projectile fragmentation with parti- 
cle beams were reported by Schimmerling, Vosburgh, and 
Todd (1971) at the former Princeton-Penn Accelerator. 
Silicon detectors provided a reliable and convenient approach 
to the measurement of fragmentation cross sections and were 
used by several groups at Fawrence Berkeley National labo- 
ratory including the work of Schimmerling et al. (1989) of 
Ne beams in a water column, Webber, Kish, and Schrier 
(1990) on carbon and CH 2 targets for cosmic-ray studies, 
and by Cummings et al. (1990) for a variety of targets of 
interest for shielding analysis. More recently a large number 
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FIG. 19. Schematic diagram showing the relative sizes, depths, 
and positions of the target and detectors used in experimental setup 
of the Lawrence Berkeley Laboratory group. The position of the 
detectors and target materials and thicknesses are modified depen- 
dent of the projectile charge and energy and the target material. 
From Zeitlin et al., 2010. 


of measurements have been made with a similar approach by 
the Lawrence Berkeley Laboratory group with the experi- 
mental setup shown in Fig. 19. Upstream of the target are 
detectors that ensure that events are initiated by a single 
incident ion. The T1 and T2 detectors are 330 mm thick 
with active areas of 300 mm 2 and are used to trigger the 
experiment. A pair of 1 mm thick with 1500 mm 2 active area 
position-sensitive silicon detectors (PSD1Y and IX) are 
placed downstream of the target material. A second set of 
position-sensitive detectors (PSD2Y and 2X) are placed 
30 cm farther downstream. Each PSD generates two 
position-dependent signals and a signal proportional to the 
total charge liberated in the detector. Data analysis examines 
the energy loss in T1 and T2 and the position-sensitive 
detectors. Based on the energy loss in the various detectors, 
a series of cuts is made to ensure the primary initiated the 
fragmentation event and to reject events where fragments 
undergo a second reaction in one of the detectors. The 
number of events of a given charge is denoted N(Z), and 
the probability of the projectile interaction in the absence of 
the target by P 0 (Z). 

The charge changing cross section is given by 
a = -Aln[P(primary)] ( 

cc pdN a 

where N a is Avogadro’s number, p is the target density, and A 
is the target mass number. For describing projectile fragmen- 
tation, total inclusive cross sections for fragmentation of the 
more abundant GCR nuclei are needed at several energies and 
for an array of targets of interest for spacecraft shielding and 
for transport in tissues. The cross section for each fragment is 
corrected for the background from the detector system 
(Zeitlin et al., 1997). 

For heavy ion fragmentation, a minimum data set includes 
elemental distributions for Z F > 2; however, isotopic cross 
sections are more stringent tests of cross-section models than 
the elemental cross sections. For neutron, hydrogen, and 
helium fragments, cross section differentials in energy and 


Rev. Mod. Phys., Vol. 83, No. 4, October-December 2011 


Marco Durante and Francis A. Cucinotta: Physical basis of radiation protection in . . . 


1265 


TABLE VI. Representative experimental data on fragmentation of Z = 10 to 28 projectiles on 
elemental targets. 


Projectile 

Energy 

(GeV /nucleon) 

Targets 

Z F range 

Reference 

20 Ne 

0.6 

H, C, Al, Cu, Sn, Pb 

3-9 

Zeitlin et al. (2001) 

20 Ne 


H, C 


Webber, Kish, and Schrier (1990) 

24 Mg 

3.65 

C, Al, Cu, Ag, Pb 

6-11 

Sampsonidis et al. (1995) 

24 Mg 

0.6 

H, C 

6-11 

Webber, Kish, and Schrier (1990) 

28 Si 

14.5 

H, C, Al, Cu, Ag, Pb 

6-13 

Brechtmann and Heinrich (1988) 

28 Si 

0.45 

H, C, Al, Cu, Ag, Pb 

6-13 

Flesch et al. (2001) 

28 Si 

0.6 

H, C, Al, Cu, Ag, Pb 

5-13 

Zeitlin et al. (2008) 

32 S 

3.65 

C, Al, Cu, Ag, Pb 

7-15 

Sampsonidis et al. (1995) 

32 S 

0.7 

H, C, Al, Cu, Ag, Pb 

6-15 

Brechtmann and Heinrich (1988) 

32 S 

1.2 

Al, Pb 

6-15 

Brechtmann and Heinrich (1988) 

35 C1 

0.65, 1.0 

H, C, Al, Cu, Ag, Pb 

5-16 

Zeitlin et al. (2008) 

36 Ar 

0.4 

H, C, Al, Cu, Ag, Pb 

7-17 

Iancu, Flesch, and Heinrich (2005) 

40 Ar 

1.65 

C, KC1 

9-19 

Tull (1990) 

40 Ar 

0.4 

H, C, Al, Cu, Ag, Pb 

5-17 

Iancu, Flesch, and Heinrich (2005) 

40 Ar 

0.4, 0.65 

H, C, Al, Cu, Ag, Pb 

5-17 

Zeitlin et al. (2008) 

48 Ca 

0.21 

Be 


Westfall et al. (1979) 

48 Ti 

1.0 

H, C, Al, Cu, Ag, Pb 

5-21 

Zeitlin et al. (2008) 

56p e 

1.09 

H, C 

12-25 

Webber, Kish, and Schrier (1990) 

56p e 

1.55 

H, C, Al, Cu, Pb 

12-25 

Cummings et al. (1990) 

56p e 

1.05 

H, C, Al, Cu, Pb 

12-25 

Zeitlin et al. (1997) 

56p e 

0.66 

H, C, Al, Cu, Ag 

6-25 

Flesch et al. (1999) 

56p e 

1.65 

H, C, Al, Cu, Ag 

6-25 

Flesch et al. (1999) 


angle are required. Table VI shows a recent survey of existing 
measurements for elemental distributions for fragments pro- 
duced from GCR projectiles with charges from 10 to 26 
showing the target mass projectile energy available. Cross 
sections for the C and O nuclei are more numerous and have 
considered isotopic distributions rather than elemental ones 
(Olsen et al., 1983). Data sets for projectiles can be divided 
roughly into three types of studies, with the first including 
more than one beam energy over the full range of target 
masses important for spacecraft and planetary habitats. A 
second series of experiments made complete measurements 
at one beam energy over the full range of target masses with 
supplemental data at other energies. Finally, other data sets 
consider measurements limited in target mass and energy 
spread. There are a reasonable number of projectile fragmen- 
tation data sets now available with the above noted excep- 
tions; however, the projectile energies of the data are lacking 
with more data needed in the 0.1 to 0.4 GeV / amu region and 
above 1 GeV / amu. Also cross-section data and multiplicities 
extending to lower fragment charge are needed. The isotopic 
composition of the GCR (Webber, Kish, and Schrier, 1990; 
Cucinotta, Kim, and Ren, 2006) needs to be considered with 
the full grid of isotopes rather than reduced ones in order to 
eliminate unnecessary error. The Fe fragmentation measure- 
ments of Zeitlin et al. (1997) resolved discrepancies in older 
data sets. Figures 20 and 21 compare measurements to several 
nuclear reaction models. The odd-even effect is not described 
by the nuclear fragmentation model version 2 (NUCFRG2) 
and is under-represented by the phits code. The QMSFRG 
model is in good agreement with many data sets (Cucinotta, 
Wilson et al., 2006) and provides a good description of the 
odd-even effects, which originates largely from the ablation 
step. 

Light particle production (n, p, d, 3 H, 3 He, 4 He, and 
mesons) in the GCR is largely from primary and secondary 


protons and neutron induced reactions. Current transport 
codes predict 10%-30% of light particles are produced by 
helium and HZE nuclei induced reactions depending on 
material type. However, these estimates are based on sparse 
cross-section databases. A large number of measurements 
exist for p and n production for p and n induced reactions 
[see, e.g., Alard et al. (1987) and review by ICRU (2000)]. 
Most radiation transport codes have assumed only an evapo- 
ration component for d, t, h, and a target fragments and 
neglecting the importance of a fast knockout component of 
these ions for GCR transport. 

These ions contribute substantially to the buildup effect in 
shielding due to their ranges which scale with the proton 
range as p:d:t:h:a to 1 :2:3:3/4: 1 and provide up to a 25% 


61 0 MeV/nudeon 40 Ar + C 



FIG. 20. Fragmentation cross sections for 40 Ar on 12 C at 
610 MeV/nucleon. Comparisons to several fragmentation models 
are shown. From Zeitlin et al., 2008. 
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40 Ar+ 12 C (0.6 GeV/u) 56 Fe+ 12 C (0.6 GeV/u) 




FIG. 21 (color). Comparisons of the QMSFRG model (Cucinotta, Wilson et al., 2006) to experiment for isotopic distribution of fragments 
for (left panel) 40 Ar on 12 C at 0.6 GeV/nucleon and (right panel) 56 Fe on 12 C interactions at 0.6 GeV/nucleon. Experimental data from 
Webber, Kish, and Schrier, 1990. 


increase in the dose equivalent for aluminum shielding over 
the evaporation contribution alone (Cucinotta, Townsend 
et al., 1996). Direct knockouts of high-energy tritons are 
especially important because of their large ranges in shield- 
ing. Knockout components will be sensitive to nuclear struc- 
ture effects such as shell structure and clustering in the 
nuclear ground state. For example, a knockouts contribute 
for 12 C and lfi O projectile or targets, d knockout for 14 N, and t 
for 27 Al. The limited existing measurements of these cross 
sections suggest that double-differential cross-section mea- 
surements at several proton energies on a wide range of 
targets are needed. Light particle production from composite 
projectiles provides an interesting effect where a substantial 
number of particles are produced with a higher velocity than 
that of the projectiles due to internal Fermi motion of the 
projectile. The effects of this process on GCR transport have 
not been studied. These data could be supplemented with 
double-differential cross sections for light particle produc- 
tion from several of the more abundant GCR nuclei (e.g., He 
and O). 

D. Validation of transport codes 

There are three approaches to assessing uncertainties in 
transport models describing exposures to sensitive tissue sites 
behind spacecraft shielding from space radiation: 

(1) Comparison of ground-based measurements for de- 
fined beams of thin and thick targets for different 
material compositions and amounts. 

(2) Intercomparison of radiation transport codes using 
matched configurations and environments. 

(3) Comparison of transport codes to spaceflight 
measurements. 

The NASA Space Radiation Laboratory (NSRL) made 
extensive measurements for a variety of FIZE nuclei of the 
Bragg ionization curve in polyethylene or aluminum shield- 
ing. A recently developed Monte Carlo based transport code, 
the GCR event-based risk model (germcode), utilizes the 


QMSFRG model of the nuclear interaction database, and the 
atomic energy loss subroutines from hzetrn to describe the 
NSRL beam line for radiobiology applications (Cucinotta, 
Plante et al., 2011). Figure 22 shows a comparison of the 
NSRL measurements to the germcode for 28 Si, 37 C1, 48 Ti, 
and 56 Fe nuclei. Excellent agreement between the model and 
measurements is seen at all depths, including past the Bragg 
peak or primary-ion range where only secondary radiation 
contributes. In Fig. 23 we show comparisons of the 
germcode to measurements by Zeitlin et al. (2008) of 
elemental distributions of secondary fragments at two depths 
in polyethylene shielding for 1 GeV/u Fe beams. Agreement 
between theory and measurements is typically within ±20%. 
Energy loss in silicon detectors has been measured behind 
graphite epoxy. Data were compared to the greentrn code, 
which has many important overlaps with the HZETRN code. 
The results in Fig. 24 show good agreement between the 
code and measurement (Walker et al., 2005). 

Intercomparisons of the hzetrn code to different 
Monte Carlo codes (hetc-heds and fluka) have been per- 
formed for well-defined shielding configurations and identi- 
cal source energy spectra for GCR (Heinbockel et al., 2011a) 
and SPEs (Heinbockel et al, 2011b). These results, as 
summarized in Fig. 25, show good agreement for neutron 
and deuterons from GCR proton spectra behind 20 cm of 
water shielding. Other comparisons on organ doses did not 
consider all of the contributions from target fragments that 
introduce minor differences and some inconsistencies 
in how the different codes evaluate the dose are apparent 
(Heinbockel et al., 2011a, 2011b). 

These comparisons show a reasonable overall agreement 
between transport code predictions when compared with 
identical source spectra and shielding configuration. 
However, in some cases different codes give very different 
predictions, especially for the production of light ions, where 
predictions can differ as much as an order of magnitude. A 
possible minor discrepancy is the mesons, electrons, and 
y rays that may contribute 10% to 15% of the absorbed 
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FIG. 22 (color). Comparisons of the GERMCODE with QMSFRG cross sections to NSRL measurements for the depth dose in polyethylene 
for nearly monoenergetic 56 Fe (0.59 GeV/u), 48 Ti (0.98 GeV/u), 37 Cl (0.5 GeV/u), and 28 Si (0.403 GeV/u) nuclei. From Cucinotta, Plante 
et al, 2011. 



FIG. 23 (color). Comparison of the GERMCODE with the QMSFRG 
nuclear cross-section model to thick target data from Zeitlin et al. 
(2008) for fragmentation of 1 GeV/nucleon 56 Fe beam at two 
depths of polyethylene. 


dose from GCR behind ISS levels of shielding (Aghara et al., 
2009); however, they provide a much smaller contribution to 
the dose equivalent, because of their small quality factor 
compared to HZE nuclei or stopping protons and He nuclei. 
These processes have not been fully integrated into all ver- 
sions of the various existing codes and should be considered 
for deep shielding predictions (> 50 g/cm 2 ). Agreement 
between the NASA codes with phantom torso measurements 
on the space shuttle and ISS are within 20% (Cucinotta et al., 
2008), and the comparison of energy spectra discussed above 
suggests a similar agreement would be obtained with other 
transport codes. 

Energy spectra for light particles have been measured with 
particle hodoscopes and shown good agreement with the 
hzetrn code as shown in Fig. 26 (Badhwar et al., 1995) 
and for neutrons in Fig. 27. Measurements of heavy ion 
spectra require large detector areas and few have been 
made in human rate vehicles because of the mass require- 
ments. Other spaceflight measurements of LET or micro- 
dosimetry measurements of lineal energy spectra using 
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Energy Loss (MeV) 


FIG. 24. Comparisons of the GREENTRN code to measurements for 
the summed energy loss in the silicon detector from fragments of a 
1053 MeV/amu 56 Fe beam behind 10g/cm 2 graphite epoxy 
(51%/49%). From Walker et al., 2005. 


E. Track-structure models 

Track-structure models provide a description of the posi- 
tion of excitations and ionization of target molecules from the 
passage of ions through a medium. These descriptions are 
needed in theoretical models of biological responses, for 
understanding the extrapolation of limited radiobiology data 
to other radiation qualities, and for describing the response 
of radiation detectors to the variety of ions in space 
(Katz, Scharma, and Homayoonfar, 1972). Application of 
Monte Carlo track-structure simulation codes has been used 
for studying the distribution and types of initial DNA damage 
including models of single-strand break, double-strand break, 
base damage, and clusters of DNA damage (Goodhead and 
Nikjoo, 1989), of oxidative species created by radiation 
(OH , H 2 0 2 , e — , etc.) (Plante and Cucinotta, 2008) or the 
description of the response of TEPCs (Nikjoo, Khvostunov, 
and Cucinotta, 2002). Analytic models have also been used 
for a range of application in particle detection and the theo- 
retical description of radiobiology data. The assumptions of 
analytical models can be validated by the more detailed MC 
simulations and are useful for space radiation applications, 
where a large number of particle types and energies place 
practical limitation on the use of MC simulations. 
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FIG. 25 (color). Intercomparison of secondary neutron spectra 
calculated with several space radiation transport codes. From 
Heinbockel et al., 2011a. 
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TEPC have also been shown to be in good agreement with 
transport codes (Shinn et al., 1999; Barber et al., 2001). 
Agreement is improved when models of detector response are 
used to make such comparisons. Corrections of measure- 
ments often involve model dependent assumptions such as 
the conversion from silicon to tissue-equivalent LET spectra 
where the conversion factor is dependent on energy and 
charge number. Also, the combined TLD plus CR-39 method 
for estimating dose equivalent (Yasuda et al., 2001) involves 
corrections on the high LET sensitivity of TLDs and low LET 
sensitivity of CR-39. TEPC comparisons to measurements 
require the conversion of energy and charge spectra into 
lineal energy spectra, which can differ substantially from 
LET spectra. 


1. Monte Carlo track simulations 

Ionization and excitation processes caused by the ions 
track and the electrons liberated by the primary ions lead to 
a stochastic process of biological events as particles pass 
through DNA, cells, or tissues (Paretzke, 1987). The ener- 
getic secondary electrons denoted as S rays can traverse many 
cell layers from the track. The so-called core and penumbra 
are loose terminology that appear in the literature; however, 
there has been no precise definition of these terms other than 
the core is close to the track and the penumbra due to the far 
reaching S rays. Older models (Chatterjee and Schaefer, 
1976) artificially assumed 50% of the energy deposition is 
in the core and 50% in the penumbra. In fact, more than 80% 
of the energy deposition from heavy ions is imparted by 
<5 rays (Cucinotta et al., 1998). 

The cross sections needed for track simulation codes are 
the total, total elastic, total inelastic, ionization, and excita- 
tion cross sections. Single and double-differential cross sec- 
tions are needed for production of electrons from the primary 
ion and for production from the secondary <5 rays (Plante and 
Cucinotta, 2008). For ion tracks, codes include ritrack, 
PITS, OREC, and DELTA; however, the energy and charge range 
of many of these codes are limited to energies below 
10 MeV/nucleon or require extensive CPU times to sample 
high-energy 8 rays produced by relativistic ions. The 
ritrack code generates radiation tracks up to very high 
energies (> 10 GeV/nucleon) and includes relativistic cor- 
rections to ionization and excitation cross sections, including 
bremsstrahlung contributions (Plante and Cucinotta, 2008, 
2009) as well as the production of different radiolytic species. 

The GSI code trax (Kramer and Kraft, 1994) uses the 
empirical cross-section formulas by Rudd (1997) modified 
for use at high ion energies and with empirical corrections in 
accordance with energy loss tables. With this correction, the 
microscopic calculations of trax can be expanded to predict 
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FIG. 26 (color online). Comparisons of proton (left) and deuteron (right) energy distributions from GCR on the STS-48 mission to HZETRN 
results. From Badhwar et al., 1995. 



spectra measurements. 

macroscopic quantities, such as the position of the Bragg 
peaks of heavy ions. Bragg peak positions can be reproduced 
within less than 0.7 mm, and the dose values are in good 
agreement with macroscopic calculations, trax yields the 
width of the Bragg peaks in a natural way as a result of the 
emission of 8 electrons with a spectrum of kinetic energies. 
With the inclusion of nuclear interactions, which are not 
described in the current code version, trax should be able 
to accurately reproduce the macroscopic dose deposition 
patterns starting from ab initio microscopic track- structure 
simulations (Kramer and Durante, 2010). 

There are few data at high energies to validate the accuracy 
of the electron production cross sections or effective 
charge assumptions used by these models (Rudd, 1997). 
Intercomparisons of codes for the frequency distribution in 
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FIG. 28. Projections over the X-Y plane of simulated tracks 
segments (calculated at ~ 1 0 12 s) for the following impact 
ions: 4 He 2+ (0.45 Me V am u '), 12 C 6+ (I MeVamu '), 28 Si l4t 
(90 MeVamu -1 ), and 56 Fe 26+ (1 GeV amu '). Ions are generated 
at the origin along the Y axis in liquid water at 25° C under identical 
LET conditions (~ 150 keV /im 1 ). Each dot represents a radio- 
lytic species. From Plante and Cucinotta, 2008. 


biomolecules by ions and electrons have been made and 
indicate reasonable agreement between models (Nikjoo, 
Uehara, and Brenner, 1997). In Fig. 28 we show a stochastic 
description of energy deposition in water for different ions 
with identical LET of 150 keV/yu,m. 


2. Analytic track-structure models 

A “first-order” numerical approach that relates LET, the 
basic physical parameter in conventional risk assessment, to 
track structure is to consider the relationship between the 
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radial dose distributions about the track and the LET. The 
radial dose is the energy density distribution in a cylindrical 
shell of radius f, about the ions path (Butts and Katz, 1967). 
LET is related to the radial dose by integrating the radial 
distribution over all radial distances up to the maximum 
allowable value t M : 

LET = 2ir [ U tdt[D s (t) + D exc (f)] + nucl. stopping. 

Jo 

(36) 

In Eq. (36) contributions from ionization are denoted here 
as D s , and excitations D exc are considered in the radial 
distribution. The value of t M defines the track width and is 
a function of ion velocity, corresponding to the range of 
electrons with maximum energy ejected by the passing ion. 
The track width can extend well beyond 100 /xm as the ion’s 
velocity approaches the speed of light. The effects of nuclear 
stopping in radiation action are only important for very low- 
energy ions ( < 0. 1 MeV/nucleon). In the model of Kobetich 
and Katz (1968) the primary electron spectrum from ion 
interactions with target atoms is folded with average trans- 
mission properties of electrons to obtain the spatial distribu- 
tion of electron dose as a function of radial distance from 
the ion’s path. The radial dose from ionization and 5-ray 
transport is described by 

(37) 


In Eq. (37) w is the initial electron energy, E is the residual 
energy of an electron with energy a> after traveling distance t, 
and r\ (t, co) is the transmission probability that an electron 
with starting energy cu penetrates a depth t. Equation (37) 
includes an angular distribution for the number of primary 
electrons produced from target atom i, n h with energy co and 
solid angle fl. The cross sections for electron production 
from protons are scaled to heavy ions using the effective 
charge Z*. The angular distribution has important effects on 
the radial distribution at both large and small radial dis- 
tances, and only a minor effect at intermediate values where 
a 1/f 2 behavior holds. An ansatz (Brandt and Ritchie, 1974) 
can be used for the radial dependence of the excitation term 
D cxc (t) as 


£>cxc« = C exc (A, Z, 0) 


exp(-t/2d) 


(38) 


where d = f3/2hc/(2ira) r ), c is the speed of light, fi is the 
ion velocity scaled by c, h is Plank’s constant, co r = 13 eV 
for water, and C exc is the normalization parameter. In 
Eq. (38) the radial extension of excitations is confined to 
very small distances ( < 10 nm) as characterized by the 
parameter d (Brandt and Ritchie, 1974). Characteristics of 
the two components of the radial dose are illustrated in 
Fig. 29, for two ions of LET close to 30 keV//xm (1 MeV 
H and 300 MeV/nucleon Ne ions) (Cucinotta, Nikjoo, and 
Goodhead, 1999). The radial dose for the Ne beam extends 
for many microns, while the low-energy proton beam depos- 
its all of its energy within 0. 1 /xm of the track. 



FIG. 29. Comparison of calculations of radial dose distributions to 
experiments for *H at 1 MeV (LET = 27 keV/yum) and 20 Ne at 
377 MeV/nucleon (LET = 31 keV//u,m). From Cucinotta, Nikjoo, 
and Goodhead, 1999. 


A “second-order” numerical approach to track structure, 
more closely related to the Monte Carlo simulation, is to 
model the frequency distribution of energy imparted to a 
volume of biomolecular dimensions (Goodhead, 1995; 
Cucinotta, Nikjoo, and Goodhead, 2000). For high-energy 
ions, the frequency distribution can be described using the 
ion’s impact parameter and distinguishing events where the 
ion passes through the volume (primary-ion events) and out- 
side the volume (5 -ray events), and by determining the mean 
and variance of the energy imparted to the volume including 
corrections for 5-ray escape out of the volume. The two 
components are weighted by considering the number of 
events as a function of impact parameter. The frequency 
distributions (Goodhead, 1995; Cucinotta, Nikjoo, and 
Goodhead, 2000) demonstrate that there are energy deposi- 
tion events in biomolecular targets that occur for high-LET 
radiation that are not possible with low-LET radiation, even at 
high doses (< 100 Gy) (Goodhead and Nikjoo, 1989). 
Validations of theoretical predictions of track-structure mod- 
els are severely limited at this time; however, the develop- 
ment of new experimental techniques in nanodosimetry could 
offer approaches for validation. 


IV. SHIELDING 

On Earth, radiation protection strategies are normally di- 
vided into 3 groups as follows: (1) increasing the distance 
from the source, (2) minimizing time of the exposure, and 
(3) use of radiation shielding. 

Increasing the distance to the source plays no role in space 
because the GCR are omnidirectional, and SPEs permeating 
along magnetic field lines become isotropic within a 
few hours after an event start (NCRP, 2006). Minimizing 
time of the exposure is problematic for space exploration 
because it may limit the performance of mission science 
objectives. Increasing rocket velocity to shorten times to 
planetary destinations is a long-term goal of space programs. 
The Apollo missions to the Moon lasted up to 12 days, 
exploration of near Earth objects such as asteroids may take 
up to a year, and a Mars mission can last up to 1 100 days with 
conventional rockets. New nuclear reactor powered plasma 
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propulsion systems such as the VASIMR are under develop- 
ment, which can cut Earth-to-Mars transit times below 
60 days (Chang-Diaz, 2000), but at the time exist in concept 
only (Durante and Bruno, 2010). For farther destinations, 
such as the outer planets or interstellar exploration, long times 
in space must be accommodated by radiation protection 
methods. Thus, developing spacecraft or planetary habitat 
shielding becomes the major approach for radiation protec- 
tion with biological countermeasures or genetic selection for 
radiation resistance outside the scope of this review. 

SPE shielding problems have been shown to be readily 
solved using current technologies and localized storm shelters 
inside spacecraft, which is not the case for the problem of 
GCR shielding. For both GCR and SPE, material selection 
and optimization of topology are major considerations. 
Spacecraft volumes may be constrained as well as mass 
when considering shielding augmentation. More importantly, 
the extra fuel required to launch such shielding compounds 
the mass dedicated to shielding. There is also a competition 
between shielding mass relative to other necessary resources 
or flight safety factors. Dual use shielding approaches, such as 
water, fuel, and food storage, are useful in this regard. In situ 
shielding using planetary resources is also of interest. Passive 
and active shielding approaches are both areas of ongoing 
research. 

Basic concepts in atomic and nuclear physics can be used 
to guide material selection augmented with detailed radiation 
transport code predictions. Engineering considerations on 
material strength, temperature, UV degradation, flammabil- 
ity, etc. must be considered along side of radiation protection, 
and the composite picture must be analyzed. Materials with 
the smallest mean atomic mass are usually the most efficient 
shields for both SPE and GCR. The composition of the 
radiation field changes as particles lose energy and suffer 
nuclear interactions in traversing structural materials, instru- 
ments, and the tissues of astronauts. Both the energy loss and 
the changes in particle fluence are related to the number of 
atoms per unit mass (in units such as grams) in the traversed 
material, which, in turn, is proportional to Avogadro’s num- 
ber divided by the atomic mass number A for each element 
of the material. The energy loss by ionization of a single 
component of shielding material with atomic number Z is 
proportional to the number of electrons per atom and thus 
proportional to Z/A. However, the energy lost per gram of 
material and per incident fluence (e.g., in units of particles 
per cm 2 ), the “mass stopping power,” is also inversely 
proportional to the density p (e.g., in g/cm 3 ) of the material, 
so that the energy lost by one incident particle per cm 2 per 
unit mass is proportional to (Z/ p)A. 

The number of nuclear interactions per unit mass and per 
unit incident fluence is proportional to a/ A, where cr is the 
total nuclear reaction cross section. To a first approximation, 
cr is proportional to A 2 / 3 , so that the nuclear transmission is 
proportional to 1/A 1 / 3 . The ratio of electronic stopping power 
to nuclear interaction transmission is thus proportional to 
(Z/ p)A 2 / 3 . Materials with small atomic mass have the highest 
number of electrons per nucleon (e.g., Z/A is 1 for hydrogen, 
0.5 for carbon, but 0.48 for aluminum, 0.46 for iron, and 0.40 
for lead). Light mass materials have smaller nuclei and there- 
fore more of them can fit into a given mass, so that there can 


be more nuclear interactions. Furthermore, the ratio of ion- 
ization energy loss to nuclear interactions is also dependent 
on the material density. For liquid hydrogen (p = 
0.07 g/cm 3 ), the ratio is —14, whereas for aluminum 
( p = 2.7 g/cm 3 ) the ratio is only 0.5, and for lead (p = 
11.3 g/cm 3 ) the ratio is 0.2. From these considerations, an 
electron plasma would provide the best shield from GCR. A 
shield made of liquid hydrogen, which has the highest ratio of 
electrons to nuclei per atom and produces minimal secondary 
radiation (e.g., mesons), is the second best choice. Hence, 
interest in polyethylene and hydrogen embedded nanofibers. 
The character of these interactions and the secondary nuclei 
produced through both projectile and target fragmentation is 
important. Lighter nuclei have fewer neutrons to release and 
some nuclei, e.g., carbon, can break into three helium nuclei 
without releasing any neutrons. In tissue, the release of three 
helium atoms is much more biological damaging than that of 
neutrons; however, if produced within spacecraft shielding 
materials neutrons are a higher concern because of then- 
longer ranges than slow helium particles. For very thick 
shields, lighter nuclei are also more effective in shielding 
against the built-up neutrons. For these and related reasons, 
detailed knowledge of the actual composition of the radiation 
fields (and of the biological consequences of exposure to 
them) is required to evaluate the net effect of shielding 
materials. 

A. Shielding of solar particle events 

Solar particles are constrained in energy with most event 
containing particles with kinetic energies below a few hun- 
dred MeV, thus allowing effective shielding approaches to be 
developed with passive shielding. As described above there 
have been about 400 SPEs in the space age, however, many of 
these events lead to very small organ doses. Kim et al. 
(2009a) ranked each event by the proton fluence above 
30 MeV to form a probability distribution of the likelihood 
of an event of a given size and time period using hazard 
analysis. They then considered the variability in energy spec- 
tra for the events with integral fluences > 10 8 protons/cm 2 , 
corresponding to 34 SPEs. Considering the variability of SPE 
spectra is preferred over the choice of a single or small 
number of event spectra (ASEB, 2008). Figure 30 shows 
results for the blood-forming organ doses behind 5 g/cm 2 
aluminum shields. This analysis shows that only 4 SPEs in the 
space era would have led to doses sufficient to possibly cause 
acute radiation sickness with fairly minimal shielding 
amounts. Furthermore, no events would have led to acute 
radiation death with minimal shielding due to their cumula- 
tive organ dose and dose rates (Hu et al., 2009). The results 
also show that the integral 30 MeV proton fluence is a poor 
predictor of organ doses, and a more recent analysis shows 
sufficient improvements when the 100 MeV integral proton 
fluence is used as a predictor (Kim et al., 2009b). The 
situation would be different for EVAs where doses can reach 
high levels for many events; however, spacesuit materials and 
design reduce exposures to some extent (De Angelis et al., 
2004). Using the catalog of past SPEs and probabilistic risk 
assessment approaches, the level of shielding needed is es- 
tablished to be in the few to about 20 g/cm 2 for the most 
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FIG. 30 (color online). The blood-forming organ (BFO) dose 
behind 5 g/cm 2 aluminum shields from the 34 largest SPEs in 
the space age (1955-2010) ranked by 30 MeV proton integral 
fluence, T> > 30 MeV calculated by the BRYNTRN code. The circles 
show the doses from each of the 34 events. The error bands show the 
90% confidence levels due to the variability of the proton energy 
spectra based on the 34 events. The dashed line shows the NASA 
30 day dose limit for the BFO. From Kim et al., 2009a. 

dangerous events. The wall thicknesses of spacecraft are on 
the order of 5 g/cm 2 ; however, when the internal equipment 
is considered shielding thicknesses much higher occur, and 
20 g/cm 2 of shielding is easily achieved in localized areas 
called storm shelters within space vehicles. Optimizing mass 
through material selection and topological considerations are 
the focus for SPE shields. Shielding approaches must be 
developed for transit vehicles, planetary habitats, and EVA 
including pressurized rovers each of which has specific limi- 
tations (ASEB, 2008). Figure 31 shows the advantages of 
polyethylene or water compared to aluminum in protecting 
against SPEs where the probability of exceeding the 30-day 
limit has been calculated for increasing 100 MeV proton 
fluence. A sufficient reduction in SPE shielding requirements 
occurs if materials with high hydrogen content are utilized. 

Because SPEs last only several hours at the peak 
exposure rates, localized shielding approaches including 
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FIG. 31 (color). Calculations of the probability to exceed the 30- 
day dose limits for the blood-forming organ of 250 mSv, for 


EVA shielding, portable shielding, and using crew sleep 
quarters as storm shelters are considered to be sufficient for 
SPE protection. Mass requirements on the order of 1000 kg 
are needed to provide storm shelter protection for a crew of 
three when optimal material selection and topologies within 
vehicles or habitats are considered. For EVA shielding nearby 
pressurized rovers is advantageous; however, spacesuit de- 
sign and shielding blankets or coats made of hydrogenous 
materials to cover vital organs is advantageous. 


B. Shielding of galactic cosmic radiation 

Figure 32 shows the contributions from various charge 
groups to the effective doses for a Mars mission near solar 
minimum predicted by the hzetrn transport code (Cucinotta, 
Hu et al., 2010). Charge group contributions to the effective 
doses for the interplanetary and surface segments are shown 
along with the total for a 30-month mission to Mars. Note the 
neutron effective dose from Z = 1 and 2 particles is counted 
in the total values for these groups and not in the neutron 
category, which shows only Z > 2 contributions from neu- 
trons. In Fig. 33 the energy spectra for different charge groups 
are shown. The heavy ions attenuate slowly with increasing 
shielding; however, slow protons and helium nuclei of rela- 
tively high-LET values build up in tissue along the path of the 
higher-energy protons, neutrons, and heavy ions making up 
about one-third of the tissue dose equivalent. The LET spec- 
tra for increasing amounts of water and liquid hydrogen 
shielding are shown in Fig. 34 (Wilson et al., 1995). 
Clearly LH 2 provides a significant reduction in high-LET 
components compared to other materials. The predictions 
of transport codes can be tested by heavy ion accelerator 
measurements. Zeitlin et al. (2008) considered a large num- 
ber of materials exposed to 1.087 GeV/nucleon Fe nuclei 
using the detector system uses to measure fragmentation 
cross sections described in Fig. 19. The measurements (see 
Fig. 35) verify the improved efficiency of the hydrogeneous 
materials. A complication arises due to the many energy and 
charge groups in space. In fact, Durante et al. (2005) showed 



s 

Radiation type 


FIG. 32 (color). Charge group contributions to the effective doses 
for Mars mission and interplanetary and surface segments. Note the 
neutron effective dose from Z = 1, 2 particles is counted in the total 
for these groups and not in the neutron category which shows only 
Z>2 contributions from neutrons. From Cucinotta, Hu et al., 
2010. 
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FIG. 33. Contributions of different GCR charge groups vs kinetic 
energy to the skin dose equivalent calculated in the hzetrn code for 
5 g/cm 2 aluminum shielding. From Wilson et al., 1995. 

that for chromosomal aberrations behind shielding some ions 
have reduced effectiveness behind shielding and others in- 
creased, dependent on where on the RBE versus LET re- 
sponse curve the particle falls. 

The importance of various GCR charge groups and ener- 
gies in shielding material selection is dependent on the 
biological response model considered. However, these con- 
siderations are clouded by the uncertainties in such knowl- 
edge. Current radiobiological response models of cancer and 
other risks are not able to determine which shielding material 
is optimal under statistical tests that consider the uncertainties 
in such a calculation (Cucinotta, Kim, and Ren, 2006). A 
more recent concern is the possibihty of nontargeted effects 
or other peculiarities of the response of humans to radiation 
(Barcellos-Hoff, Park, and Wright, 2005; Cucinotta and 
Chappell, 2010) which lead to a deviation from a linear 
response model; the value of radiation shielding and the 
importance of mission length are diminished if biological 
responses are sublinear with increasing doses and possibly 
increased if threshold doses occur. 

Slowing down incident GCR particles using materials with 
a preponderance of energy loss due to ionization and a 
minimum probability of nuclear interactions is not always 
an optimal strategy. Nuclei such as carbon and oxygen, 



FIG. 35. Dose reductions for various materials placed in a 
1.087 GeV/nucleon 56 Fe beam. From Zeitlin et al., 2008. 


incident at high energy, have low LET, well beneath the 
peak value of RBE. When they lose energy in a shielding 
material (without suffering nuclear interactions), their LET 
increases. As a consequence, their RBE also increases instead 
of decreasing, so that they become more hazardous, rather 
than less hazardous. On the other hand, the LET of incident, 
high-energy heavier nuclei, such as iron, is close to the 
150 keV//r,m corresponding to the peak RBE. Losing energy 
and slowing down further increases their LET beyond the 
peak, yet they may become no more hazardous despite their 
higher LET due to the role of cell killing. Nuclear interactions 
that change a penetrating GCR nucleus into lighter nuclei, 
e.g., nuclear interactions that fragment Si into C and He 
(a particles), result in particles of lower RBE, the desired 
outcome. On the other hand, fragmentation of oxygen or 
carbon atoms in tissue into a particles produces particles 
near the maximum RBE of energies of just a few MeV. The 
fragmentation of high-LET iron into lower-LET chromium or 
sihcon fragments changes the contribution to the radiation 
field from less hazardous particles, beyond the RBE peak, to 
more hazardous particles with an LET at or before the peak. A 
complicated picture emerges with many balances between 
increasing and decreasing energy deposition, and secondary 
projectile and targetlike nuclear fragments. 

1. Shielding performance index 

Aluminum structures have been the main stay of space- 
craft. Also clear is the fact that hydrogen is certainly the 
optimal material for radiation protection. Materials with 





FIG. 34. Differential LET spectra from the HZETRN code for increasing amounts of liquid hydrogen, water, and aluminum shielding. From 
Wilson et al., 1995. 
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atomic constituents heavier than aluminum such as steel or 
lead produce too many neutrons to be of interest for space 
radiation protection. With these ideas in mind, Wilson et al. 
(1993) devised a shielding performance index to bound ma- 
terial effectiveness relative to aluminum when considering 
different biological response models. Shield performance of a 
material m is defined relative to the standard aluminum 
technology by the shield performance index, which is the 
ratio of a risk assessment variable in aluminum to that of 
material m. Figure 36 shows the contributions from different 
GCR charge groups at solar minimum comparing the con- 
ventional systems using radiation quality factors to radio- 
biological models of cell killing, neoplastic transformation, 
and mouse tumors. Using these results, Fig. 37 shows the 
performance index relative to increasing aluminum shielding 
amounts. Performance indexes from 2 to more than 10 are 



FIG. 36 (color). Several different measures of radiation quality 
calculations of different GCR charge groups to overall risk measures 
for 5 g/cm 2 of aluminum and no tissue shielding are shown. From 
NASA, 1998. 



FIG. 37. Range of maximum performance index relative to liquid 
hydrogen relative to aluminum using various biological response 
models including Harderian gland tumors in mice (HG), trans- 
formation (T), and dose equivalent (H). 


possible with typical amounts of spacecraft shielding, 
although the calculations used only a minimal amount of 
tissue shielding and thus likely overestimate the perfor- 
mances. These results suggest that new high hydrogen con- 
tent materials could provide significant benefits compared to 
aluminum, especially if amounts from 10 to 30 g/cm 2 could 
be afforded in vehicle designs. 


2. Uncertainty analysis 

An uncertainty analysis using Monte Carlo propagation of 
errors from different factors that enter into a shielding calcu- 
lation of cancer risks was developed by Cucinotta, Manual 
et al. (2001) and Cucinotta, Kim, and Ren (2006). 
Uncertainty factors considered are epidemiology data, dose- 
rate effects, radiation quality factors, and space physics. 
Probability distribution functions (PDF) representing the 
range of possible values for each factor are devised. The 
cancer risk projections for space missions are found by 
folding predictions of the tissue-weighted LET spectra be- 
hind spacecraft shielding dF/dL with the radiation cancer 
mortality rate to form a rate for a trial J: 

r dF 

mj(E, a E , a) = m u {E, a E , a) / dL — LQ lmi _j(L)x L -j 


(not shown are the uncertainty quantiles associated with the 
low-LET mortality rate, mu). The result of Eq. (39) is then 
used with the formula for the risk of exposure-induced death 
(REID) that considers a double-detriment life table for radia- 
tion cancer risk on top of the background of other causes of 
death in the U.S. population. After accumulating sufficient 
trials (~ 10 5 ), the results for the REID estimates are binned 
and the median values and confidence intervals found. A \ 2 
test for determining if PDFs for two distinct shielding con- 
figurations or materials are significantly different was devel- 
oped (Cucinotta, Kim, and Ren, 2006). They denoted the 
calculated PDFs for a REID of R t for two configurations or 
materials as Pi(R,) and p 1 (R i ), respectively. Each p{R t ) 
follows a Poisson distribution with variance -Jp{R t ). The x 2 
test for characterizing the dispersion between the two distri- 
butions is then 


M S. (40) 

" pjiRn) + pl(Rn) 

Once x 2 i s determined, the probability P(n, x 2 ) that the 
two distributions are drawn from the same parent distribution 
is calculated. If x 2 is sufficiently large such that Pin, x 2 ) is 
less than about 20% is an indication that we can conclude that 
the two distributions lead to distinct cancer risks from GCR 
and/or SPEs, with the material with the lowest mean and 
upper 95% Cl (confidence interval) values preferred for 
radiation protection. However, the opposite result indicates 
either that the materials are approximately the same or that 
the uncertainties in risk models prevent us from concluding 
that the configuration or material is superior for radiation 
protection properties. Cucinotta, Kim, and Ren (2006) eval- 
uated x 2 f° r the LET-dependent parts of the uncertainties 
(quality factors and physics) separately, since only these 
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FIG. 38 (color). PDFs for fatal cancer risk probability for 40-year- 
old males on a 600-day Mars swing-by mission with 20 g/cm 2 
shielding of aluminum, polyethylene, or liquid hydrogen. Effective 
doses and point estimates or 95% Cl for risk of exposure-induced 
death (REID) are shown in the box. From Cucinotta, Kim, and Ren, 
2006. 


contributions explicitly depend on the modification of radia- 
tion fields by shielding. 

Figure 38 shows PDFs for the GCR as solar minimum for 
aluminum, polyethylene, and liquid hydrogen shielding with 


areal density of 20 g/cm 2 . Table VII shows results for the x 2 
statistic comparing polyethylene and hydrogen shielding ef- 
fectiveness to aluminum. The August 1972 SPE is considered 
for these comparisons. As shown by Fig. 30, large variations 
will occur in risk estimates for specific events; however, the 
August 1972 SPE provides a good model for a large event. 
For calculations that ignore the low-LET risk coefficient 
uncertainties, near solar maximum with a large SPE, signifi- 
cant results (P < 0.01) are found for liquid hydrogen, and for 
the GCR at solar minimum (with P < 0. 15). We conclude that 
because of the modest differences between polyethylene and 
aluminum as GCR absorbers, and the large radiobiological 
uncertainties in cancer risk projection models, the benefits of 
polyethylene compared to aluminum shielding for GCR can- 
not be proven at this time. A key factor in these results is that 
although aluminum is a greater producer of secondary radia- 
tion compared to polyethylene, tissue shielding “equalizes” 
much of these differences. 

3. Novel and in situ shielding materials 

Polyethylene has been identified as a useful structural 
polymer for spacecraft shielding with various fabrication 
strategies developed for damage tolerant stiff structure and 
inflatable vehicles (Wilson et al., 1991). An important chal- 
lenge is how to form high-density polyethylene fiber with 
polyethylene matrix and bond the resulting composite face 
sheets to form polyethylene foams. Approaches considered 


TABLE VII. From Cucinotta, Kim, and Ren (2006). x 2 test at solar minimum for 20 g/cm 2 shields for 40-yr males on a Mars swing-by 
mission. P(n, x 2 ) is the probability materials cannot be determined to be significantly different (n = 500). Values in bold for Pin, x 2 ) < 0.2 
indicate a significant improvement over aluminum. 


Test material 

E (Sv) 

REID (%) 

95% Cl 

X 2 /n 

P(n, x 2 ) 

Aluminum 

0.87 

3.2 

All uncertainties 
[1.0, 10.5] 



Polyethylene 

0.78 

2.9 

[0.94, 9.2] 

0.05 

>0.99 

Hydrogen 

0.43 

1.6 

[0.52, 5.1] 

0.63 

>0.99 

Aluminum 

0.87 

3.2 

LET-dependent uncertainties 
[1.9, 8.7] 



Polyethylene 

0.78 

2.9 

[1.8, 7.5] 

0.08 

>0.99 

Hydrogen 

0.43 

1.7 

[1.0, 4.2] 

1.10 

<0.15 

Test material 

r solar maximum for 5 g/cm 2 
E (Sv) 

shields. 

REID (%) 

95% Cl 

X 2 l n 

P{n, X 2 ) 

Aluminum 

1.21 

4.4 

All uncertainties 
[1.5, 13.1] 



Polyethylene 

0.94 

3.5 

[1.2, 10.8] 

0.14 

>0.99 

Hydrogen 

0.52 

2.1 

[0.60, 6.4] 

0.81 

>0.99 

Aluminum 

1.21 

4.4 

LET-dependent uncertainties 
[3.0. 11.0] 



Polyethylene 

0.94 

3.5 

[2.3, 8.8] 

0.32 

>0.99 

Hydrogen 

0.52 

2.1 

[1.2, 5.2] 

1.38 

<0.001 

Test material 

r solar maximum for 20 g/cm' 
E (Sv) 

1 shields. 

REID (%) 

95% Cl 

X 2 ln 

Pin, x 2 ) 

Aluminum 

0.54 

2.0 

All uncertainties 
[0.60, 6.8] 



Polyethylene 

0.45 

1.7 

[0.52, 5.6] 

0.08 

>0.99 

Hydrogen 

0.24 

0.9 

[0.27, 2.9] 

0.77 

>0.99 

Aluminum 

0.54 

2.0 

LET-dependent uncertainties 
[1.2, 5.6] 



Polyethylene 

0.45 

1.7 

[1.0, 4.6] 

0.15 

>0.99 

Hydrogen 

0.24 

0.9 

[0.52, 2.4] 

1.26 

<0.005 
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are plastic thermosets, aliphatic systems, and e-beam curing. 
Aliphatic systems lose strength at both low and high tem- 
peratures to be expected under space conditions. Wilson, 
Cucinotta, and Schimmerling (2004) discussed the usage of 
aliphatic and aromatic hybrid polymers as an improvement to 
a purely aliphatic system. Carbon nanotubes with high hy- 
drogen content offer a distinct approach and are under 
investigation for improved structural layouts. 

The high costs of launch mass make in situ shielding on the 
lunar or Mars surface of interest. Distinct approaches may be 
required since a higher relative fraction of astronaut expo- 
sures on the Mars surface are due to neutrons compared to the 
lunar surface due to the Mars atmosphere both depleting 



Energy (MeV) 


FIG. 39 (color). Neutron doses on Mars. Variations in charged 
particle doses are caused by the atmospheric height in different 
geometric locations, whose maps were measured by the MOLA 
altimeter on the Mars Global Surveyor spacecraft, launched in 1996. 
The upper panel shows the seasonal variation of the integral neutron 
fluence above 10 MeV due to ice buildup reducing regolith back- 
scattered neutrons. The bottom plot shows the forward and back- 
ward neutron fluences for different Mars soil compositions. From 
Wilson, Clowdsley et al., 2004. 


heavy ions and adding neutrons. The exposure to astronauts 
will depend on which region of Mars a habitat occupies 
(Wilson, Clowdsley et al., 2004), because large variations 
in secondary neutrons can occur due to the soil composition 
as well as seasonal variations due to the presence of protec- 
tive ice (water or permafrost) as shown in Fig. 39. Landing 
sites will likely be chosen by science requirements; however, 
they may require distinct shielding approaches dependent on 
the fractional contribution from neutrons to the total risk. 
Digging regolith to shield planetary habitats with take-along 
equipment is one approach considered. Certainly, exposures 
could be reduced substantially if several meters of regolith are 
used to cover a habitat. However, long times spent within 
habitats will likely reduce the science objectives of such 
missions. The Moon is known to contain lava tubes as well 
as craters that partially shield space radiation, and Mars has 
similar geographical features to increase habitat shielding. 
In situ hydrogen or perhaps launched boron could be used to 
create regolith shielding bricks; however, they will require 
new manufacturing approaches on planetary surfaces. 


C. Active shielding 

Active shielding involves the generation of electromag- 
netic fields to deflect space radiation. Several approaches 
have been proposed (Townsend, 2001) as follows: (1) electro- 
static shield, (2) plasma shield, (3) confined magnetic field, 
and (4) unconfined magnetic fields. 

Each approach should be considered relative to passive 
shielding in terms of effectiveness as well as power and mass 
requirements and the possibility of device failure. Also, harm- 
ful effects need to be considered including biological effects 
of strong magnetic field on humans, and the disturbing of lunar 
dust, which contains nanophase iron (ASEB, 2008) from the 
electromagnetic field. Lunar dust was reported to irritate the 
skin and respiratory systems of the Apollo astronauts. 

Developments in superconducting magnets are necessary 
to provide effective active shielding from GCR. The experi- 
ence gained during the development of the alpha magnetic 
spectrometer (AMS) superconducting magnet has been useful 
to develop ideas and techniques to be applied to radiation 
shield for exploration missions (Battiston, 2007, 2008). 
However, an anomalous heating of the cryogenic supercon- 
ducting magnet eventually forced the AMS Collaboration to 
exchange the superconducting AMS-02 magnet (Blau et al., 
2004) to the previously developed, less sensitive, permanent 
magnet. This permanent magnet has been delivered to ISS by 
the Space Shuttle Endeavor on 16 May 2011, and is now 
operating in its search for antimatter. During the development 
of the cryogenic superconducting magnet (2 coils of Ni-Ti 
operating at 1.8 K, generating a field of 0.87 T in the center), 
the AMS Collaboration proposed conceptual designs of mag- 
netic configuration optimized for a travel to Mars or for a 
Moon base, with toroidal configurations able to provide a 
zero-magnetic field inside the spacecraft. However, the 
experience of the AMS Collaboration also shows that the 
cryogenic superconducting magnet technique is critical for 
spaceflight and presently not yet reliable. In the future, 
cryogen-free superconductors could provide excellent solu- 
tions, pending technological improvements in this field. 
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As an example magnetic lenses built from superconducting 
magnets have been described by Spillantini et al. (2007) for 
protection from directional SPE-type radiation. An angular 
cone ±10° could protect a several cubic meter volume 
against protons with energies up to about 200 MeV using a 
so-called D-shaped toroid field with a devise mass of about 
1100 kg, which is in the range of the mass required by a 
passive shielding storm shelter to protect the same volume. 
However, the directionality of the field is a disadvantage 
compared to passive shielding. 

Recently, renewed interest in electrostatic shielding has led 
to several possible designs. The model described by Tripathi, 
Wilson, and Younquist (2008) illustrates many of the prob- 
lems with such approaches. Tripathi, Wilson, and Younquist 
(2008) described a lunar base electrostatic shield made of 12 
spheres: 6 of 20 m radius located 160 m from the habitat 
along several axes, and a 6 of 10 m radius located 50 m along 
50 m away from the habitat. Each sphere carries a +300 MV 
charge. The resulting configuration was shown to reduce the 
GCR dose equivalent at the surface of the skin by about two- 
thirds. A smaller reduction would be found at the deep seated 
organ due to the buildup of target fragments from the residual 
high-energy GCR nuclei. The mass of shields approached that 
of the habitat and therefore comes at a very high launch cost, 
perhaps of interest for a permanent moon base. These possi- 
bilities will require very large power requirements, and nu- 
clear powered reactors are a useful option. 

V. SPACE RADIATION EFFECTS 
A. Biological effects 

Most of the uncertainty on space radiation risk is associ- 
ated with poor knowledge of the biological effects of cosmic 
rays. In particular, gaps in knowledge are mostly related to 
the following: 

• relative biological effectiveness factors of energetic 
heavy ions for late effects, both cancer and noncancer; 

• dose and dose-rate reduction effectiveness factors; 

• errors in human data including statistical, dosimetry, 
and transfer between populations in application to space 
radiation risks; 

• effects of the exposure to the mixed high- and low-LET 
space radiation field; 

• shape of the dose-response curve at low doses for 
charged particles; 

• and interaction of radiation damage with other space 
environment stressors (particularly microgravity). 

The main biological effects associated with exposure to 
cosmic radiation are as follows: 

• carcinogenesis, 

• late degenerative tissue effects, 

• acute effects, and 

• hereditary effects. 

Cancer is currently dominating risk estimates, and as dis- 
cussed the dose limits are most constraints by cancer mortal- 
ity risk. However, noncancer effects are becoming an 
increasing source of concern. They can again be divided into 

• acute and late damage to the CNS, 
cataract formation, 


• cardiovascular diseases including coronary heart disease 
and stroke. 

• digestive and respiratory diseases, 

• and accelerated senescence leading to endocrine and 
immune system dysfunction. 

Acute effects are expected only for very intense SPE 
hitting an unprotected crew. The prodromal risks such as 
nausea, vomiting, and fatigue are the most likely (Hu 
et al., 2009), and acute mortality is very improbable. 
Hereditary effects are generally assumed to be a factor of 
10 or more lower than somatic effects in radiation protection, 
although there is missing evidence for heavy ions in this field. 
The age at flight and the small number of astronauts tend to 
reduce the importance of hereditary and reproductive effects. 

Interestingly, so far the cataract is the only cosmic 
radiation-induced effect actually observed in astronauts 
(Cucinotta, Manual et al., 2001; Chylack et al., 2009). 
Recent epidemiological evidence does not support the exis- 
tence of dose threshold for radiation-induced lens opacifica- 
tion, although the shape of the dose-response curve is still 
unclear (Blakely et al., 2010). However, animal studies show 
that the RBE of heavy ions for cataract induction is as high as 
50 at doses below 100 mGy (Brenner et al., 1991), thus 
explaining the high effectiveness for accelerated caratacto- 
genesis in space. 

The RBE represents a major source of uncertainty, and 
indeed the problem of radiation quality is central in radiation 
protection on Earth as well. The shape of the dose-response 
curve at low doses is a key problem, and the importance of 
nontargeted effects or adaptive response mechanisms is still 
unclear both for radiation protection on Earth (Mullenders 
et al., 2009) and in space (Cucinotta and Durante, 2006; 
Cucinotta and Chappell, 2010). Conventional risk assessment 
approaches have been limited to using averaged quantities 
such as dose or dose equivalent. The stochastic models of 
nuclear reactions and track structure described in Sec. Ill will 
be needed to fully utilize the results of radiobiology experi- 
ments and to develop new stochastic based disease models to 
replace the oversimplifications embodied in the conventional 
RBE approach. 

The large number of health risks of concern, entailing 
many tissues and diseases, complicates pursuing approaches 
such as genetic selection, because it is highly unlikely that 
one person could be resistant to so many diseases. 
Approaches to biological countermeasures are also chal- 
lenged by the multitude of health risks as well as the types 
of radiation in space and the chronic exposure of up to 3 years 
for a Mars mission. These factors limit the applicability of 
radio-protectors developed for terrestrial exposures, which 
are usually for acute exposures to low-LET radiation. 

Large ground-based radiobiology research experimental 
programs are currently ongoing to reduce uncertainty on 
biological effects of heavy ions, in the USA (at the NASA 
Space Radiation Laboratory in the Brookhaven National 
Laboratory, NY), Europe (at GSI in Darmstadt, Germany), 
and Japan (at HIMAC in Chiba). Results from these experi- 
ments have been reviewed (NCRP, 2000; Durante, 2002; 
ESA, 2006; NCRP, 2006; Cucinotta and Durante, 2006; 
Durante and Cucinotta, 2008; NASA, 2009; Maalouf, 
Durante, and Foray, 2011), and the interested reader will 
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find all of the most updated biological information in those 
papers. 

B. Effects on microelectronics 

The effects of cosmic radiation on microelectronics are 
also a major concern for the safety in spaceflight. In fact, the 
charge generated by the ionization of energetic charged 
particles in spacecraft electronics can upset microcircuit 
functions and eventually lead to catastrophic failures. 

HZE damage is caused by electron-hole pairs that are 
generated in silicon dioxide and other insulators (McNulty 
et al., 1996; Holmes-Siedle and Adams, 2002). Proton- and 
neutron-induced single-event effects are generally caused by 
nuclear reaction fragments produced in the semiconductor, 
typically spallation nucleons and Si recoils (Tang, 2005). 
The damage in the integrated circuit is the result of trapping 
of excess charges at or near the interface region between the 
oxide (or other insulator) and semiconductor. The effect de- 
pends on the particle charge and energy, and on how excess 
electron-hole pairs are transported within the oxide. When low 
electric fields are present in the oxide, recombination reduces 
the pairs transported to critical interface regions (charge yield). 
If the charge yield is low, then the net effect of the deposited 
dose is reduced. If the charge yield is high enough, it can 
induce destructive effects (single-event latchup, single-event 
gate rupture, single-event burnout, single-event dielectric 
rupture, and single-event snapback) or nondestructive effects 
(single-event upset, single-event transient, single-event func- 
tional interrupt, and single-event hard errors). 

Single-event upset (SEU), or “soft errors,” are the most 
common space radiation-induced effects, and they occur also 
in our terrestrial computers. It is a change of state (bit flip) 
affecting bipolar and metal-oxide-semiconductor elements. 
As noted, SEU are nondestructive, i.e., the bit flip can be 
corrected by rewriting the affected element, but they can still 
represent a threat to spacecraft electronics. Microelectronic 
failures or multiple-bit error rates generated serious problems 
in several missions, including Cassini and the Hubble tele- 
scope. Since they represent a major issue in the design of any 
kind of space mission, these effects have been reviewed 
several times [see, e.g.Jniewski (2010), and the reader can 
consult the literature for specialized information]. What is 
particularly interesting from the physical point of view, and in 
a way similar to the problem of biological effects, is the 
dependence of SEU from the energy and charge. 

Before electronic components can be sent into space, they 
have to be qualified on Earth. These tests are usually done 
using y rays (for total dose effects) or at low-energy accel- 
erator facilities (for SEU). The vulnerability of a device to 
SEU is normally defined by a threshold LET (minimum LET 
required to produce SEU) and the saturated cross section 
(maximum upset rate for the device, when no increase is 
observed by further increasing the LET). But, as for the 
biological effect, LET alone may not be the only parameter 
characterizing the effect. Below the LET threshold, nuclear 
reactions can still produce SEU, and high-energy ions are 
more effective than low-energy ions (Dodd et al., 2007). On 
the other hand, recent studies using particles at energies 
around 1 GeV/nucleon suggest that above the LET threshold 


and before saturation, the cross section is higher for low- 
energy ions compared to relativistic HZE at the same LET 
(Connell etal., 1996; Metzger et al., 2011), a track-structure 
effect (see Sec. III.E) also observed for inactivation of viruses 
(Butts and Katz, 1967) and mammalian cells (Katz et al., 
1971; Cucinotta, Wilson et al., 1996; Elsasser et al., 2010) 
by heavy ions. Tests at high-energy accelerators are also 
preferable because the particles can cross the whole device 
and allow testing the device at different tilt angles (Pellish 
et al., 2010; Metzger et al., 2011). High-energy spallation 
sources can also be used for neutron-induced SEU testing. 


VI. CONCLUSIONS 

Along with biomedical research, efforts in nuclear, parti- 
cle, and radiation physics are essential to the success of future 
exploration missions to Mars and other destinations. There 
are many areas where physics takes a lead role in risk 
assessment for astronauts, including the translation from 
radiobiology and biophysics into radiation protection meth- 
ods and spacecraft design. Radiation physics has shown that 
shielding should be successful for SPE protection and how to 
optimize approaches in order to reduce costs and improve 
safety factors. Physics has also taken the lead in determining 
how poorly shielding approaches against GCR work, based 
on the detailed understanding of nuclear reactions and sec- 
ondary radiation that occur from the passage of protons and 
HZE nuclei in all materials including tissue. Analysis also 
indicates that active shielding using magnetic fields and 
electrostatics is severely challenged for GCR, yet continues 
to be of interest because of the potentially high payoff. 

Existing radiation transport codes agree with measurements 
of organ dose equivalents within 20%; however, further im- 
provements in transport codes are of interest. Future work will 
continue to improve knowledge of interaction cross sections 
and their usage in space radiation transport codes. Many cross- 
section data sets have been determined; however, there are still 
significant gaps in double-differential cross sections for light 
particle production, cross sections for new materials of inter- 
est, and knowledge of exclusive cross sections that is needed to 
support stochastic based risk models that supplant the current 
deterministic dose based quantities. The description of energy 
deposition in biomolecules due to nuclear processes or S rays 
as mapped to tissue and organ structures is needed to support 
new risk assessment approaches as well. In parallel, robotic 
precursor missions to the surface of Mars and other near Earth 
objects sample the environment with particle detectors and 
dosimetry prior to the arrival of future explores. Transport 
codes and energy deposition models play essential roles in 
the understanding of these new data sets. 
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